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CHAPTER I 


GENERAL INTRODUCTION 


1,1 INTRODUCTION 

"The term air pollution may be defined as any undesirable 
addition to our environment in the form of dust# particulate 
matter# fume# gas# mist# smoke, smog etc#, for duration and in 
such quantities that it becomes harmful to human, animal and 
plant life, to our living, to our enjoyment and to our property" 
Air pollutants (compounds of sulfur, carbon, nitrogen, etc#) 
which are emitted from industrial and other man made sources 
affect the ecosystems in general and its specific components 
in particular in several ways leading to harmful consequences* 

It affects our health, our living, our property, our priceless 
monuments etc* It also harms animal and plant species and 
damages the ecosystem creating various kinds of problems in the 
form of diseases, epidemics, climate changes, etc# 

It is ^therefore# absolutely essential to study the 
dispersion of pollutants, emitted by various industries and man- 
made projects, rec[uirad for environmental assessment programme 
of polluted regions on one hand and to determine its impact on 
both living and nonliving species on the other# 

In view of this, in this thesis the following type of 


oroblemsare investiaated 
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(i) Unsteady state dispersion of reactive air pollutant 
forming secondary pollutant in the atmosphere with 
various kinds of removal mechanisms including deposition 
on the ground 

(ii) Steady state dispersion of reactive pollutant with 
variable wind velocity and diffusivities from point 
and line sources by considering deposition on the 
groxind. 

(iii) Steady and unsteady state dispersion of air pollutants 
in presence of green belt* 

In the following wa present an outline of the 
relevant literature so that the work done in the thesis can be 
seen in its proper perspective# 

1*2 OUTLINE OF RELEVANT LITERZvTURE 

In the study of pollutant dispersion# we are generally 
interested in finding out the concontra tion of a pollutant, 
after being emitted into atmosphere# at a specified location 
and time, which depends upon the n\mbar of sources, their 
locations, emission rates, etc# The dispersal process is 
governed by the process of molecular diffusion and convection 
and dependsupon meteorological factors such as wind, temperature 
inversion, topography of terrain, removal mechanisms such as 
dry deposition, rain out/washout, etc# (Calder, 1961/ Pasquill, 
1962/ Turner, 1970/ Carpenter et al*,i971/ Heines and Peters, 
1973a, i973b/ Whaley, 1974* Morgenstern at al* 1975/ Seinfeld, 



3 


1975; Shum et al»^ 1975; Erniak^ 1977; Peterson and Seinfeld,1977/ 
Dobbin S/ X97^f Fisher and Macqueen/ 1981; Novotny and Chesters / 
198ljr Wark and Warner, 1981; Karamchandani and Peters, 1933) • 

In particular, the solution of steady state three dimensional 
diffusion equation with constant diffusivities and wind velocity 
has been obtained by incorporating settling velocity, dry 
deposition on the ground (Ermak, 1977) ♦ Fisher and Macqueen 
(1981 ) have solved cross wind integrated atmospheric diffusion 
equation under the inversion condition* 

In the last few decades, the increased number of studies 
have been carried out to understand the scavenging of contaminant 
by rain (Postma, 1970; Hales, 1972, 1973, 1932a, 1932b; Slinn, 
1974; Chu and Seinfeld, 1975; Scriven and Fisher, 1975a, 1975b; 
Peters, 1976; Maul, 1973; Barrie, 1973; Fisher, 1932; Kumar, 

1985, 1936)* In particular^ Postma (1970) and Hales (1972) have 
presented some fundamentals of the precipitation scavenging of 
gases by rain* Hales (l972) has also called attention to the 
possibility that a gas plume could become redistributed in the 
atmosphere under the action of reversible washout* Slinn (1974) 
has analytically considered the redistribution of gas pl\xme 
caused by reversible washout- Scriven and Fisher (i975a, 1975b) 
have discussed the long range transport of air borne material 
and its removal by deposition and washout* An excellent review 
of the removal processes and models of long range transport of 
air pollutant has been presented by Fisher (1983)# A theoretiral 
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analysis of contribution to rainwater sulfate concentration by 
precipitation scavenging of gaseous S 02 and sulfate (So^ ) 
aerosols has been presented by peters (1976)* An Eulerian 
model for simulating the coupled processes of gas phase 
depletion and aqueous phase accumulation of pollutant species 
has been investigated by Kumar (1935)* Kumar (1936) has also 
extended the above Eulerian model for rain scavenging of 
pollutants by taking into account the process of absorption of 
multiple pollutant species and chemical reaction within the 
rain drops* The combined effects of dry deposition on the 
ground and various removal mechanisms such as rainout/wa shout, 
chemical reaction etc# on the dispersion of pollutants by 
considering them as first order process have also been studied 
by many workers (Astarita et al#/(1979; Alam and Seinfeld, 1931)* 
In particular, Alam and Seinfeld (1981 ) have discussed the 
dispersion of sulfur dioxide and sulfate from continuous point 
source and obtained analytical solution of steady state three 
dimensional diffusion equation by considering their removal due 
to rainout/wa shout process and dry deposition on the ground# 

Peters and Richards (1977) have developed a procedure to 
study the steady state dispersion of pollutants frcra a point 
source with simultaneous reversible reaction when the ground 
surface acts as a reflecting plane# The validity of the 
procedure necessitates the chemical reaction to be rapid enough 
relative to the transport so that local chemical equilibrium is 
achieved at the receptor site# Shukla et al#, (1932) have 
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discussed the first order reverisible absorption of a pollutant 
from an area source in the presence of fog layer* 

The dispersion of air pollutant frciii line source has 
also bean studied (Walters, 1969, Dilley and Yen, 1971/ 

Liu and Seinfeld, 1975, Sharma and Myrup; 1975, Chock, 1978, 1932; 
Hassid, 1933). In particular, Walter (1969) has studied the 
dispersion by considering uniform air velocity and taking eddy 
diffusivities as a linear fianction of height* Chock (1973) has 
proposed a simple linear source model to describe the downwind 
dispersion of pollutants near the roadways* 

It is noted here that the sources of pollutant could 
be time dependent (for example, explosion of nuclear device 
such as reactor accident at Chernobyl /U.S.S.R., leakage of 
toxic gases from a plant as happened at Bhopal, India) and 
little attention has been paid to find exact solution of 
atmospheric diffusion equation in such cases particularly 
with removal mechanisms such as chemical reaction and deposition 
on the ground (Cleary et al* 1974* Llewelyn, 1983)* 

In view of this, in chapters II to V, the 
dispersion of reactive air pollutant from time dependent source 
by incorporating first order chemical kinetics and forming 
secondary pollutant has been discussed* In these chapters, wind 
velocity and diffusivities are taken as constant* 

In general, the wind velocity and diffusivities are 
functions of space coordinates and time. In the case when 
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horizontal homogenity is asstamed^ the wind velocity and diffusion 
CO efficients may be taken as functions of height only* For 
general fimctional forms of these variables no exact solution 
of three dimensional diffusion equation with chemical kinetics 
and other removal mechanisms is known* In few cases, however, 
when the wind velocity and diffusion coefficients are functions 
of height, only, the analytical solution of steady state two 
dimensional diffusion equation for non reactive air pollutant 
has been obtained when the ground surface acts as a reflecting 
plane (Smith, 1957* Demuth, 1978; Robson, 1933) » Heines and 
Paters (1974) have presented the solution of atmospheric diffusion 
equation (when wind speed is constant and diffusivities are 
function of downwind distance) by incorporating dry deposition 
on the ground in the model# As stated earlier, xjnder horizontal 
homogenity, the meteorological variables are functions of 
height only# Since a general function can be approximated 
by a set of step functions, one way to study this kind of 
problems is to divide the inversion height into several layers 
such that the diffusion coefficients and wind velocity are 
assumed constants but different in different layers* The 
diffusion equation with these constant values is then solved 
in each layer with appropriate matching conditions at the 
Interfaces of these layers* Using this concept, in Chapter VI, 
the steady state dispersion of reactive pollutant emitted from 
point and line sources is studied under inversion condition 
incorporating dry deposition on the groiind by dividing the 
inversion layer into two, and three layers# 
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At this jiinctiore^ it may be pointed out that the 
pollutant in the atmosphere can be removed by introducing 
suitable artificial removal mechanisms to protect a given 
region* The mechanisms could be in the form of green belt 
plantation/ spraying of liquid drops including water at a 
required place and time* It is well toown that certain 
plants have the capability of removing pollutants in the 
atmosphere by absorption, deposition etc* and this property 
can be used for protecting a given region (Petit et al*/l976i' 
Bache, 19793/ 1979bf Smith/ 1981 jr Kosher et al*/1932? Slinn, 
1982)* In particular, the transport and deposition of particulate 
matters within plant and vegetable canopies have been studied* 
Petit et al* (1976) presented some results concerning 
characteristics of air flow within and above a forest by 
calculating S 02 fluxes at the top of the canopies* Bache (1979 )ha 
also developed an analysis and suggested a modified form of 
diffusion equation to study particulate transport within and 
above foliage canopy* Slinn (1981) gave a theoretical frame 
worh to predict particle deposition due to vegetation by 
considering wind velocity profile* A review of atmospheric 
deposition regarding plant assimilation of gases and particles 
has been presented by Smith (1981) and Kosher et al* (1982)* 

A mathematical model for aerosol depletion and deposition on 
forests has also been suggested by considering a modified form 
of convective diffusion equation (Wiman, 1985)* This model 
considers the interaction between forest structure and open 
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field by considering forest aerodynamics and aerosol 
cha rac t eri sties* 

Frcan the above studies it may be speculated that if a 
suitable green belt in the form of forest is provided aroxmd 
the pollutant source or some distance away from it and close 
to the place (habitat, historical monximents, etc*) to be 
protected, then it is possible to reduce the concentration of 
pollutant by this green belt and thus protecting the region 
under consideration# In this direction Kapoor and Gupta (1934) 
studied the attenuation of an inert gas by green belt under the 
steady state condition# 

Keeping in view of the above, the steady and unsteady 
state dispersion of a reactive air pollutant in the presence 
of green belt is studied in Chapter VII and VIII# 

The removal mechanisms due to absorption of air pollutant 
by fog or rain droplets present in the atmosphere are very 
complex# In Chapter IX, a modified form of atmospheric diffusion 
equation is suggested in the form of an integro-partial differen- 
tial equation* The effect of instantaneous and delayed removal 
mechanisms on the dispersion of air pollutant frem an elevated 
time dependent point source is investigated# 

1#3 SUMMARY 


This thesis consists of nine chapters and deals with the 
steady and unsteady state dispersion of air pollutants fron time 
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dependent sources/ forming secondary pollutants/ xandergoing 
various removal mechanisms such as rainout/wa shout/ dry 
deposition etc» 

In Chapter I, an outline of relevant literature is 
presented so that the research work done in the thesis can be 
seen in its proper perspective* 

In Chapter II/ the xonsteady state dispersion of a 
pollutant from a time dependant elevated point source under- 
going first order irreversible chemical reaction and forming 
a secondary pollutant is discussed by assuming that these 
pollutants are removed by various removal processes such 
as rainout/wa shout/ etc* The following forms of time dependent 
sources are considered : (i) instantaneous emission (ii) constant 
emission (iii) step function type emission. 

In the case of instantaneous emission, it is shown that 
the concentration of both the species (primary and secondary) 
decreases as time increases at a particular location and the 
point of maximum in the concentration distance profiles moves 
away from the source as time passes* However/ in the case of 
constant emission, the concentrations of both the species 
increase as time increases and tend to their steady state values* 
In the case of step function type emission (i*e* emission of 
pollutant occurs only upto time i't is noted that for t £ t^ 

the concentrations of both the species increase and for t > t^ 
the concentrations decrease as time increases* 
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By using the method of images^ solutions are also 
obtained when the ground is acting as a reflecting plane*. It 
is noted here that the solution thus obtained generalises the 
Gaussian puff and plume models for reactive pollutant in an 
isotropic atmosphere* 

In Chapter III, the unsteady state dispersion of air 
pollutant emitted from an elevated time dependent point source 
undergoing first order reversible and irreversible reactions 
and forming a secondary pollutant is studied* It has been 
shown that the concentration of the primary pollutant is greater 
than the secondary pollutant when the forward reaction rate is 
less than the backward reaction rate* However, the reverse is 
the case when forward reaction rate is greater than the backward 
reaction rate* When the forward reaction rate is equal to 
backward reaction rate, the concentrations of both the species 
tend to their equilibrium values which are equal* 

In Chapter IV, the unsteady state dispersion of ' air 
pollutant from a time dependent point source located at height 
h above the ground and forming a secondary pollutant is 
investigated by taking into account the effect of settling 
velocity, dry and wet depositionc (without inversion condition)- 
The analytical solution of three dimensional diffusion equation 
is obtained by considering wet deposition as the first order 
process* It is shown that as the settling velocity increases, 
the concentration at a point located below the source height 
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increases but it decreases if the point is located above the 
source height* It is also noted that the effect of removal 
process is to decrease the concentration of both primary and 
secondary species* 

In Chapter V, the same problem is studied under the 
inversion condition* The comparison has been made between the 
two cases (with and without inversion) and it is found that the 
concentration of pollutant at a point with inversion condition 
is greater than its corresponding value in the case of without 
inversion condition. A particular case of the analysis presented 
in this chapter is applied to MIC leakage from Union Carbide 
Factory in Bhopal, India* This analysis provides useful informa- 
tion regarding impact assessment of the dispersion of MIC in 
Bhopal * 

In previous chapters, it has been assumed that the wind 
speed and diffusivities are constant in the all directions but, 
in general, the wind velocity and diffusivities are functions 
of space coordinates and time* In the case when horizontal 
homogenity is assumed, the wind velocity and diffusion coeffi- 
cients may be taken as fumctionsof height only* Since a general 
function can be approximated by a set of step functions, one 
way to study such problems is to divide the inversion height 
into several layers such that the wind velocity and diffusion 
coefficients are assxiraed constants but different in each 
layers, The diffusion equation with these constant values is then 
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solved in each layer with appropriate matching conditions at 
interfaces of these layers# Using this concept, in Chapter Vi, 
the steady state dispersion of reactive pollutant emitted from 
point and line sources is studied under inversion condition by 
taking into account the dry deposition on the ground, dividing 
the inversion height into two, and three layers# The effect 
of stepwise variation of the wind velocity (by taking diffusion 
coefficients same in each layer) is to decrease the concentration 
of pollutant near the groiand as number of layers increases* The 
effect of stepwise variation of the diffusion coefficients (by 
taking wind velocity same in each layer) is to increase the 
concentration of pollutant near the ground as ntmber of layers 
increases . 

In Chapter VII, the steady state dispersion of air 
pollutant from a source in the presence of green belt is studied* 
It is observed that the effect of green belt is to reduce the 
concentration of pollutant but the reduction is more if the 
source height is less than the average height of green belt# 

In Chapter VIII, the unsteady state dispersion of reactive 
air pollutant in the presence of green belt is investigated# 

Here the whole region in the wind direction is divided into two 
regions, the green belt being in second region lying away from 
the source* It is shown that the concentration of pollutant 
decreases due to the presence of green belt* 
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In Chapter IX, a general mathematical model for 
atmospheric dispersion, in the form of an intsgro-partial 
differential equation is suggested where the integral term 
represents the removal mechanism due to absorption of air 
pollutant in fog/rain droplets present in the atmosphere# This 
model is applied to study the effects of instantaneous and 
delayed removal mechanisms on the dispersion of air pollutant 
from an elevated time dependent point source# It is shown that 
the concentration of the pollutant along the central line 
decreases due to the removal process and its magnitude in the 
case of delayed removal mechanism may be greater than the case 
of an instantaneous removal mechanism# 
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CHAPTER II 


UNSTEADY STATE DISPERSION OF AIR POLLUTANT FROM AN 
ELEVATED TIME DEPENDENT POINT SOURCE FORMING 
SECONDARY POLLUTANT 

2 INTRODUCTION 

The dispersion of air pollutant from a point source 
has been investigated by many workers using Gaussian plume 
and other models under various conditions (Smith, 1957* 

Lamb and Seinfeld, i973?Heines and Peters 1973a, i973b; 
Dobbins, 1979/ Karamchandani and Peters, 1983/ Llewelyn, 

1983)* In general, dispersal process is affected by meteoro- 
logical conditions in the environment and various removal 
process such as chemical reaction, depositions, rainout/washout 
etc* (Slinn, 1974/ Seinfeld, 1975; Novotny and Chesters, 1981) 
For example, during dispersion SO 2 is converted to sulfate 
(SO^ ) aerosol in the atmosphere and both SO 2 and SO^ are 
removed by wet and dry depositions. Effects of these removal 
mechanisms on dispersal of air pollutants have been studied 
by considering them as first order processes (Scriven and 
Fisher, 1975; Nordland 1975; Sander and Seinfeld, 1976; 

McMohan et al., 1976; Prahm et al*, 1976; Sheih, 1977; 

Calvert et al», 1978; Slinn, 1980; Alam and Seinfeld, 1981; 
Hale, 1932; ia«wtflyp,1933) . 
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In particular^ Alam and Seinfeld (1981) have studied 
the steady state dispersion of sulfur dioxide frcm a 
continuous point source forming sulfate as a secondary 
pollutant by taking into account removal mechanisms such 
as deposition* In general, however, the sources of 
pollutant are time dependent and this aspect must be 
considered in the dispersal process* In such a case \ansteady 
state diffusion equation with chemical reaction must be 
solved under suitable initial and boundary conditions depending 
upon the nature of the soiirce (point source, line source, and 
area source! 

Therefore, in this chapter, the dispersion 
of a pollutant fron a time dependent elevated point source 
undergoing first order chemical reaction and forming a 
secondary pollutant has been investigated by assxaming that 
these pollutants are removed by various processes such as 
rain out/washout etc# The following forms of flux 
at the point source have been considered in the subsequent 
analysis : 

(1 ) source with instantaneous flux 

(2) source with constant flux 

(3) source with step function type flux# 

By using method of images, the solutions have also been 
obtained when the source is located at (0,0,hg) and ground 
is acting as a reflecting plane* 



21 


2 #2 MATHEMATICAL FORMULATION 


Consider the dispersion of a reactive air pollutant P, 
from a time dependent elevated point source in the environment 
undergoing first order chemical reaction and forming a secondary 
pollutant Q both of which are renoved by mechanisms such as 
rainout/washout, etc* It is assumed that conversion of P into 
Q and their removal are first order processes with rate constants 
kg and kp as shown in the following scheme. 



S 


where S denotes the removal of both species fron the atmosphere* 


Taking the point source as the origin of the coordinate 
system (x/y/z), where x is in the direction of wind, the 
concentration of gaseous pollutant P is governed by the following 
atmospheric diffusion equation. 


3C . „ 9 C 






a^c 


ay'^ az*^ 


) - (k + k ) c 
9 


where u is mean wind velocity assumed to be constant, 
constant diffusion coefficient and t is the time* 


( 2 * 1 ) 

D is the 


The initial and boxindary conditions for (2*l) can be 
prescribed as 

c(s,t) =0, t*0 for all s = v/^^ty^+z^ > O 


( 2 - 2 ) 
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C(s,t) = 0 ass-*<» t>0 (2*3) 

“471 s^D = W(t) as s •* 0 t > 0* (2-4) 

as - 

The last boundary condition implies that the point source 
has time dependent flux W(t) which is considered to be of the 
following forms 

(1) W(t) = 6(t) 

(2) W(t) = W (constant) 

c (2.5) 

(3) W{t) * 0 < t < t 

= 0 t > t^ 

where §(•) is Dirac-delta function and is a constant* 

Similarly^ the equation for the concentration C of the 

ir 

secondary species Q can be written in a coupled form (by 
noting the material balance of the species)^ as 


2 2 

a c a c 


-r- + U r — 

at ax 


D(— + — jE + — ^) + ]cc - k C 


p p 


(2*6) 


If there is no direct emission of the secondary pollutant 
then the boxmdary conditions for equation (2*6) can be 
prescribed as 


Cp(s*t) =0 at t ss 0 

C (s#t) = 0 ass'»‘"*t>0 

p « — 



(2#7) 
(2 .S) 

(2 »9) 
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2 *3 METHOD OF SOLUTION' 

Following Astarita et al^ (1979) and Alam and Seinfeld 
(1981), equations (2 •! ) and (2*6) can be written in the 

compact form as 



' c 

ktk 0 

L 


g 


c 

-k k 


p 

P 

1 

2 2 


f c 1 


( 2 # 10 ) 


where L 


aE “ IS 


+ D(“«y + — +■ — is an operator# 

3x ay az 


To uncouple this system^ let us define a matrix 


(2 >11 ) 


L-" >J 

which can be written as a' = rmr , where 


1 

0 

k+k„ 

0 1 

1 

0 ' 

k„“k-k„ 

, M = 

M i 

g 

0 

II 

-k 

k -k-k 

1 

L p g 

J ' 

J 

p g 



If we define a new concentration B such that 

Fcl fc 1 


= R 


B ! 
L J 


(2#i2) 


the equation (2 4*10 ) can be written in uncoupled form as 

LC ^ (k+kg) C = 0 (2 >.13) 

LB - kp B =0 (2 #14) 

and from equation (2*12) we get the following relations 

B = n - (2.15) 


B = a 


k C 


to find once B and C' are determined^ 


Keeping in view equations (2 #2) - (2 #4) and (2*7)- 

(2 *9), the boxindary conditions to determine B can then be written a 
B(s,t) = 0 at t = O (2.16) 
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B(s#t) = 0 ass-»‘»^t>0 


(2.17) 


~4jis^D II 


-k 


w(t) 


as 


S 0# t > 0 


(2.18) 


To obtain the solution of (2.13) and (2*14) subjected to 
initial and boundary conditions we proceed as follows. Taking 
Laplace transform of equation (2*13) with respect to t and 
using initial condition (2*2)/ we have 


- me 


ac _ , a^G . 

u -T-r = Dt — >5- t — 5- + 

^ ax*^ 37^ az 

mm 

where C is Laplace transform of C, m 
variable. 


(2 *19) 


p+ktk and p is Laplace 
y 


The boundary conditions becane 

C(s,p) = 0 as s -• <» ( 2 * 20 ) 

-4irs^D = W(p) as s - 0 (2. 21) 

where W(p) is Laplace transform of w(t) • 


Assuming the solution of equation (2*19) to be of the 

form 

c = expC^^) f(s/p) 

and using aquations (2 *20 - 2*21) we get 

r =c pvr. r 22 - 4. S) 1/2 -1 (■> 79) 

^ 47^ e3<P L 2D ^ ® J • 12 .22 ; 

4D 

Taking inverse Laplace transform of equation (2.22)^ 

the solution of equation (2*13) satisfying the initial and 
boundary conditions can be written as : 
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C(x,Y,z*t) = W(p) ej<p [pt^-|§ - (-^ + g)^ s] dp 

(2-23) 


where y is a real positive niamber such that all the singularities 
of the integrand lie on the left hand side of the line Re(p) = y 
in Bromwich contoxir# see fig* (2*1)« Similarly, the solution 
of (2*14), satisfying initial and boundary conditions (2*16 • 

2*17) can be written as 


JLoo ^ 

B(x,y,z,t) W(p)e:q> [pt+^ s]dp 


■" 1,1 


y— ioo 


4D* 


(2*24) 


-k 


where A = = p -h kp 


Using Bromwich contour, the above integrals can be 
evaluated and we get the following expressions for concentration 
distributions of both the species for various formsof W(t) ■ 


(i) Instantaneous flxrc l*e* W(t) =s 6('t) 


C(x,y, z, t) 


w 


O rUx ^ 


(47rDt) 


ITT m ^ 


4Dt-' 


(2*25) 


AW ■ „ -.2 2 

B(x,y,z,t) = —— ^- 372 - exp[|g + l^)t - |^ ] (2.26) 


(2) Constant flux W(t) = 
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exp (^ 2 ) 

C(x/y, z,t) - {exp [-( 5 )^'^^ s] 


B(x/y,2^t) = 


1 f i ^~ut . /U-Kl /2 , , 

‘ ~ J - e sxn(-fp-) s du} 
71 b ^ iJ 


4jiDs 


(2 *27) 


u-b. 


171 -ut . ^l^l /2 

- * / re sxn(-7r-^) 


71 V, u 
‘>1 


D 


s du} 


u 


u‘ 


(2 .28) 


where b = + k -t* kg # + kp. 

It is noted that expressions (2*27) and (2*28) for 
C and B can be obtained by integrating with respect to t the 
corresponding distributions given by equations (2.25) - (2*26) 
between 0 and t # with =s W^, 

( 3 ) Flux as step type function i.e# W(t) = W ^ < b 


s 0 


t > 0 


C(x,y, z* t) 




47lDs 


{e ° [i-HCt-t^)] 


^ ow ^ 4 /-n lit 

- I / ^ sin(^^)^'^^s [l-e ° H(t-t^)] du} 

b 

(2.29) 

A expO -4)1/2 


B(x,y^Z/t) * 


47IDs 


{e ° [l-H(t-t^)l 


•1 '* -( ^■,4. u--b< i /o 'Ut 

I / i e sin(^)l/23 j-^.^ o H(t-t„)] du) 


(2.30) 
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where is Heaviside fiinction defined by Carslaw and 

Jaeger (1941) as 


H(t-t„) =0 t < 

= F+6 

o 

= 1 t > t^+e 

where e can be made as small as we please* 


The concentration C (x/y,z,t) is obtained from C(x,y/Z,t) 
and B(x#y,z,t) for each of the above cases using following 
relation 


3c C(x,y/Z,t) 

C tx,y,z^t) * B{x*y,z,t) + ^ — 

P g 

where B and C are given by equations (2 #25 — 2 #30)# 


(2.32) 


To obtain the corresponding solutions when the source 
is located at (0#0/h ) above a reflecting plane surface (at z=o) 

3 Cl 

such that ~ = 0 at z = O# we use the method of images for 

a 2 

superposition of solutions as diffusion equation is linear^ 
(Carslaw and Jaeger 1959^ Dobbins, 1979 Chapter, 9). The 
solution of various cases are, thus, given as follows : 


( 1 ) 

C^i(x#y,z,t) 


Instantaneous flux 
W 


(47lDt) 


-m 


exp [ bt - ) ] £ e 


(z-h^)2 


'4Dt 


(zth^)' 


+ e 


4Dt 


(2.33) 
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AW 


. 2.2 


( z - h ^)- 


“ T’ ^ ^ exp [ 5 "^ “ b.j ) ] fe 


(47iDt) 


37J 


4Dt 


2D 


4Dt 


(zfh^)2 

s 


+ e 


4Dt 


} 


(2*^34) 


(2) Constant flux 


C 2 j^ (x^y,z,t) = 


W _ exp (^) ^ “ 


c ™ r * 2 D ' 


4 jrDsi 


{e 


1. I ; i e-“t 


" b “ 


w . 


— T:=frs le J - e sin(-=r-; 


rb\l/2 


47rDs, 


71 b U 


S 2 < iu } 


(2.35) 


AW exp(22) 

B 2 i ( x , y ., z , t > jnBs ' " ' — ‘® 

1 


_ i 7 t e""*" sin(^)^''^s 

71 / U ® D ^ ®1 


du} 


AW ^ exp (^) 


. : c ,— 


47rDSn 


u-b. 


-h S h sin(-^)^/2 


1 


(2-36) 


(3) Step function type flxix 


C 33 ^( x , y ^ z , t ) 


w., exp(2S) -(^) 


b^l/2 


C D i [^.H(t-t )] 


47rDs- 


{e 


- ^ ^ e*"^^sin(^=^)^'^^Sj^ [l-e ° H(t-t^) 3 


ut 


dul 


W exD(~) 

° ^ 2°- te ^ 2 [i.H(t-t )] 


47lDSr 
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sin [i—e ^ du} (2*37) 

b 


B3j|^(x,y,z,t) = 


A W^(exp(|g) ^ -(-1)^/2 g 


47IDs- 


£e ° ^ [l-H(t-t^)] 


“ I / ^ H(t-t^)] du} 

^ A exp(g§) [l-H(t-t^)] 


47rDs^ 


“ H u sin( "g -— [l-e ° H(t~t^)] du} 

*"i 

(2 ♦ 38 ) 

where s^ = x^+y^ + (z~hg)^i = x^+y^+(z+hg)^ • 

It is found that equations (2*37) and (2*38) can 

also be obtained by following relations 

^ t-t 

t o 

C 3 - (x#y/ 2 ,t) = / C..(x#y, z,t')dt'-H(t-t ) J (x,y, z# t' )dt' 


ut 


0 


t-t 


(2 .39) 


Boi(x,y,z,t) = I B-. (x,y^ 2 ,t')dt'-H(t-t ) / B^l^(x,y,z,t dt 
JX 0 o 

(2 .40) 


and = W^* 
o c 

Here the concentration C can be obtained from equation 

ir 

(2*15) by sxibstituting and (i = 1,2,3) in place of C 

and B respectively. 
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For the steady state constant flux case# equations 

(2«35 - 2*36) reduce to following forms 

W exp(^g) “W ®2 

C2^(x,y,z,t) = ^ ^ + ] (2*41) 


AW exp(^) - ^ D ^ h g ^ D ^ ®2 

B,, (x,y.2,t) = - t ^ r 1 


(2 .42 ) 

y^+(z-h^)^ y^+(2+h^)^ (k+k )D 

Taking 5 << 1^ n << 1 — << 1 


kpD 

and T-j— << 1, then the linearized form of solution can be 
u 

obtained as 


»T « J-V 

C,-t (x^y, 2 ^t) = exp [ -(- 7 f- 2 )x - 




U(z~h_) _ U(2+h^) 

{exp [ + exp [ 4 dF~^^ 

. W k„x 


A W kpx 

B,,(x*y,z,t) = exp [ - -fi- - k 


U 4Dx 


U(z-h„)' 


U(z+h^)' 


{exp [ jf ^ — ]+expC 7 nl" ' - ] ) (2.44) 


To compare present model with Gaussian plume model 

(which is conmonly used), we make further transformations 

according to the relation, 

02 = ^ 

U 


and we get 



3i 


W k+k 2 

C2^(x,y^z,t) = — Sj— exp [-(-rjj--2)x - ^ 

2710 u 2 a 


(z-h 



fe * e 


(z+hg)2 
2? 


} 


A W kpx 2 

B 2 ^(x,y,z,t) = exp [- -^ - ^] X 

2 710 U 20 


-(z-ii —(z+h 

s s 

2 2-= — 

-20^ ^ « 20^ , 

t e + e } 


(2 *45) 


(2*46) 


Putting k = kg =s lip = 0, it is straight forward to 
reduce itie SB expressions to standard form for a continuous 
Gaussian plime model in the isotropic media* It may be 
observed that the solution for Gaussian puff model to non 
reactive pollutant has been extended to a reactive pollutant 
forming secondary pollutant in an isotropic media* 


2*4 RESULTS AND DISCUSSION 


To see the effects of various parameters on the 
concentration distributions of both the species in each case, 
we use following dimensionless quantities* 


r - ^ 4- 

hi 


X 


X 

s 


- y - z r 
y = / z = , k 

s s 


h^ k 

s 

D 


k IT 


k h-^ 
D 


U = 


U h^ 


w 


c r 


''s«T 


w 




D 

O 

W h^ 
c s 


(2 .47) 
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Equations (2*33 ■* 2*34) can be written in dimensionless 
form for each case as follows (dropping bars for convenience): 

(l) Instantaneous flux 


W 2 2 

exp [~ - bt - (^^-^|-)1 ^ 


- (z"l) 

4t 


+ e 


4t 


} 


(2 wl8) 


B,,(x,y,.,t) = exp [P - b,t - 1^] - 


(47lt) 


£e 


_ (z-l)‘ 


( z+1 ) ^ 


+ e 


4t 


(2*49) 


where b and b^ (dimensionless form) are 


U 




b = X- 1, + kg , b^ = 
(ii) Constant fl\ix 


*^21 ~ ^ C^j^(x*y, z, t' ) dt^ 


(2*50) 


^21 ~ ^ B^^(x^Y/Z,t' ) dt' 

o 


(2*51) 


(iii) Step function type flxix 
t 


t-t 


^31 “ ^ Cj^j^(x#y, z,t' )dt' - H(t-t^^) / C^j^(x#y, z,t' ) dt' 


(2 .52) 
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t 

S (x,y, z, t' )dt' 


t-t 

o 

H(t-t ) / B. - (x#y#z,t' ) dt' 

o o 

(2.53) 


To see the effects of various parameters on the 
dimensionless ground level concentration, the e3q;>ressions for 
C and G given by equations (2 *48 - 2*53) are calculated, by 

i? 

choosing the following set of values of parameters ; 

= 1,0, U = 5,0, k = 0,11 # kg = 0,0, 0.03, 
k = 0*0, 0.8, y = 0»0 and depicted in figs. (2. 2 - 2 •8). 

xr 

It is observed that the concentrations C, decrease 

as downwind distance increases* It is noted that as p 

increases C decreases but increases in all cases at a 

P 

particular time and location (see fig. 2»6)i, 

When the flux is instantaneous at the source, the 
dimensionless concentrations C and C are shown in figs# (2 .2 - 

XT 

2.3) for different values of t and W^ = 1#0# It is noted from 
those figures that as time increases both C and Cp decrease 
and the point of maxima in the concentration distance profile 
of each species moves away from the source# 

For the case of constant flux, the ground level 
dimensionless concentrations C and C are shown in figs# 

ir 

(2*4 - 2*6) for different values of t# It is seen from these 

fias# that both C and C increase as time increases and would 

P 

eventually tend to their steady state values* 
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For the case of step function type flux, the ground 

level concentration is plotted in figs* (2»7 - 2*3)* These 

figs* show thet concentrations c and decrease as time 

P 

t increases for t > t^* However for t £ t^, the profiles 
for C and are similar to the case of constant flux* 

The effect of removal mechanisms on C and C can be 

ir 

seen by conparing the figs* 2*2 and 2 #3, 2*4 and 2*5, 2*1 and 
2-3 * It is observed that C and decrease as and kp 

increase respectively. 
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CHAPTER III 


UNSTEADY STATE DISPERSION OF AIR POLLUTANT FROM AN 
ELEVATED TIME DEPENDENT POINT SOURCE UNDERGOING 
FIRST ORDER IRREVERSIBLE AND REVERSIBLE PROCESSES 

3#1 INTRODUCTION 

The dispersion of air pollutant from a time dependant 
point source landergoing first order irreversible chemical 
reaction has been discussed in Chapter II* Since there are 
situations in the atmosphere where the formation of secondary 
pollutant is accompanied by a reversible reaction, this aspect 
must also be studied (Peters and Richards, 1977, Shukla et 
al * , 19 32 ) * 

Therefore^ in this chapter, the xinsteady state 
dispersion of a reactive air pollutant emitted from an elevated 
time dependent point source ^xndergoing first order reversible 
and Irreversible reactions forming secondary pollutant have 
been discussed* 

3 *2 MATHEMATICAL FORMULATION 

Consider the dispersion of a reactive gaseous pollutant 
P, anittad from a time dependent elevated point source, under- 
going first order reversible reaction forming secondary 
pollutant Q and both P and Q react irreversibily with species 
present in the atmosphere* The scheme of their reactions can 
be described as follows* 
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where k^z k^, k^ are rates of reversible and irreversible 

chemical reactions, k / k are the removal rates cf P and Q due 

g P 

to removal mechanisms present in the atmosphere such as rainout/ 
washout* 

The concentration of gaseous pollutant P is governed by 
the following diffusion equation under the assumption that the 
eddy diffusivities are unifoim in all directions (= say) 

II + U II = D(|-| + - (k+k +k ) C 4- h C (3.1) 

where U is the wind velocity taken along the x-direction and 
Cp is the concentration of secondary pollutant Q. 

Taking the source as origin of the coordinate systan, 
the initial and boundary conditions are given by 

, _ X +y +z > 0 ( 3 .2 ) 

C(szt) = 0 ass-*«>t>0 (3.3) 

-47Is^D = w(t) as s -* O t > 0. (3.4) 

ds — 

The last boundary condition implies that the point source 
has time dependent flux W(t). The following foimsof W(t) are 
considered for subsequent analysis : 
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(i) Instantaneous flux at the source 

W(t) = Wq 6(t) 

(ii) Constant flux at the source 

W(t) = W ( 3 . 5 ) 

c 

(iii) Step function type flux at the source 

W(t) = 0 < t < t 

c — o 

= 0 t > t 

o 

Under similar assumptiors the transport equation 
governing the concentration of the secondary pollutant Q 
can be written as 

3c 3C 3^c a^c a^c 

4- u = D ( 4. —-- 2 . 4- — ^)-l*kC-(k +)c +k )C (3 #6) 

3^ 3X p 33^^ r b p p 

where D is diffusion coefficient of secondary pollutant* 

Sr 

The initial and boundary conditions for C can be 

ir 

written (assxming that there is no direct emission of secondary 
pollutant) as 


C (s.t) = 0 

P ' 

0 

fl 

for s > 0 

( 3 , 

. 7 ) 

Cp(s,t) = 0 

as s -» 

<» t > 0 

( 3 , 

. 3 ) 






^Jls'^D = 0 as s 

- 0 t 

> 0 * 

( 3 . 

► 9 ) 


P 3 s 
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3.3 METHOD OP SOiUTION 

Assuming that primary and secondary pollutants are 
dispersed in the atmosphare in same manner^ in the following 
analysis we take = D* In such a case^ taking Laplace 
transform of equations (3.1) and (3*6) with respect to t., we 
get 

- 2 - 2 - 2 - 

u ^ - (k+k^+k +p)c + k^ (3.10) 

^ BY d7 ® ^ P 

30 3^0 3^0 

U = D(— + “” 2 ^ ^ — ■■y^)“’(kj^+k +k ■i*p)c tkC (3.11) 

° ax ay Bz^ P P 

where C and are the Laplace transforms of C and 

respectively, p being Laplace variable. 


The boundary conditions for C and beccme 


C( s,p) 

» 0 as s -♦ oo 

(3*12) 

-4JIS^D 

3C 

= w(p) as s - 0 

(3»13) 

Cp(s,p) 

= 0 as s -» oo 

(3-14) 

-47IS^D 

— = 0 as s 0 

as 

(3.^15) 


where W(p) is Laplace transform of W(t) « 

Assuming that the concentration distributions for 
both the species are of the following forms 

(~) 

C = e fj^(s/p) (3.16) 
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Cp = e f 2 (s,p) 

the equations (3#lo) and (3*11 ) become 


(3.17) 


) s 

( D ' ^ -tn 2 ) s 


where 



f2(s,p) 4- ^ 

k+kg+hg+p 

D / 


s f 2 (s/p) = 0 
s f^(s,p) = 0 


"^2 



D 


(3*18) 

(3.19) 


and D' = . 

Solving equations (3*18), (3*19) and using the 
corresponding boxandary conditions, the e3<pressions for C and 
Cp have been obtained as 


W(p) 

~ JjiDS 


M) ^ - ( 


p 

n 


2M 


+ 


ra'+n'_ M^i/2_ 
^ s 


4 “ 


fXR'-n* 

— 

2M 


4- M) 


Ux 



m^tn^ M)l/2 
2D D 


s 

} 


Ux _/P m'+n' _ M^i/2_ 

- W(p) k r Id "''D “2D ” D*^ ® 

S " 4iDS 2M L® 


Ux _^P ^ m'+n' 

2D ^D “2D 
e 


+ d) 
^ D^ 


Mvi/2 




(3.20) 


(3*21) 


m' 


4D 


4- k + k_ 4- k_ 


where 
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M = H-k 

Taking inverse Laplace transform of expressions (3 *20) 
“ (3#2l) by using Brcamwich contour (Carslaw and Jaeger, 1941)/ 
the following expressions for concentration distribution for 
both the species are obtained in each case as follows : 

(i) Instantaneous flux . 

W ■ b 2 

C(x,y,2,t) [55 - - 4 Dt ] 


2M 


e 3 <p [2J3 “ b't - ] } 


s_ 

4Dt 


(3-22) 


W 2 

Cp(x,Y,z,t) = ^ 




4Dt 


(3.23) 


(ii) Constant flux (Continuous source) 


w exp(^) b s 

C(x,y,z,b) = ,e 


u-b/ 


1 , 1 -ut , i^l/2 , , 

~ / - e sin (— g — ) s du} 

b' 


,b'l/2 

+ ^ I / ie-"* au)] 


2M‘ 


b' 


(3 #24) 
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Cp( X/Y/ 2, t) 


”c k r = 

47rDs 2M L 


.iT i sl„(:^)i/2s au 

,b'xl /2 

- e ^ + - / - e“^"^ sin (- V-r- -.) du ] 


71 u 


(3.25) 


(iii) Step function type flux 

p^2 ^ 

C(x.y, 2 .t) = Ljm 


(l-H(t-t^)) 


. 00 ^ ^ 

1 , i "t.t r -t « o 
— / =• e i-e 

TT n ■ ^ 


u-b' . /, 

H(t-t^) ] sin(— du} 


,b'vi/ 2 ^ 
b3 -(if) s 

^ 2M 


(l-HCt-t^) ) 


- 0° ^ *. 
i i [l-e ' 

H b. ^ 


[l-e ° H(t“t )] sin(i~-)^'^^s du} ] 


(3«26) 


Cp(x,y,2,t) 


T 7 


47 X 0 S 


(i-HCt-t^) ) 


, ^ U-b/ . /„ 

. i / i e”^"^ [l-e ° H(t-t )] sin("^- - )^'^^s du 

H U 

- e ^ (l-HCt-t^)) 


« ^ _,4- 

i , 1 — Ut 

+ — / a e I 1— e 
^ 71 u L-^ 


[l-e ^ ( t-t ) ] a'in ^ s du} 


(3.27) 
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where b' = --jS— + M , b. 


m'+n' 
■"l = "T~ 


M 


b2 = - + M,b3 = + M, 


To obtain the corresponding solution when the source 

is located at (o,0,hg) above a reflecting plane (at 2 = 0 ) 

3C 

such that = 0 at 2 = 0, we use the method of images for 
superposition of solutions since diffusion equation is linear 
(Dobbins, 1979; Carslaw and Jaeger, 1959) » 

The solution for various cases are ; 


(i) Instantaneous flux 


W„ expC^S) b, -bit 


c,(x,y,z,t) = - g . "a ... e 1 [e 


4Dt 


+ e 


°2 

4Dt 


+ ^ [e' + e" ] I 


2M 

exp(^^) 


(3.28) 


[UX) -b't 

1 j-g 


4Dt 


+ e 


L- 

4Dt j 


- e 


-b' 




% -2 
jDt ^ “jDt 3 j 


(3*29) 


(ii) Constant flux 


^2(^,7, 2, t) = 


TT «/^lvl/2 

«c ph ° s_ ~ 

“ '-2M ^s7 ■*“ « 




4xrD 


_i_ 7 i e““^ au 

7IS- 1/ U ® D ’ 

^ h 

ii” ^ u sln(-K^)^^^ 

71 S 2 U 


g ) 
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b. g ‘d =2 

-3 + 2 _ 

2M *s^, S2 

. i sln(!^)^''^s. au 

ns^ u D i 

_ ^7 i du}] (3.30) 


7IS-, ■ U - D 

b' 

TT / -.ri.>i/ 2 ^ 

^ «c®-p'r5> b rr® “ 1 , e ° 


47TD 


2 M i-‘s 


.,i_7 i au 


"=2 b' “ 


D 


- “ - a. u— b/ . /o 

.^1 i ,i„(-^)l/2s, du! 


"®2 bj “ 

f^a /2 

“^D ' 


-£■ 


_ ^"d* ®2 

4 > r 


- -i- 7 1 e-^-t ,in(2;&l)2/23^ 3^ 


b^ ^ 


D 


- 7 h sin(i^^)^^^S 2 du}] 


71 S 2 j,/ u 


(3.3i) 


(iii) Step function type flux 

rUx 


Co(x#y, z,t) 


1 

-(%i/2s -.(S.)2s 

W, .,we‘“’ *i,,‘D^=2 


4jiD 


[jI {(i-H(t-t„)){ 


o s, 


£ 5 . 


} 


, 00 ^ . ut U“bf - y„ 

J. iD-f 


* “ ^ _, 4 . , u-b/ - /« 

i- / i [l-e ^(t-t )] sin( ■■ ■ A ) ^/2 

7 is^ ,_# u >- o -’ D 2 


"^^2 b' 
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b 1 ^ D ®2 

4- ^ f ( ^ i- 

* 2M ‘ ' s Sj 


) (l-H(t-t^) ) 


i_- j- 1 [l-e ° HCt-t^)] sinC ^j^— du 


b' ^ 


. ^i_7 I 6"“*= [l-a“''° H(t-t„)] du)] 


ns2 t,, u 


Cp3(x/y,z*t) = 


Wc exp(|g) ^ 
47TD 2M 


{(• 


K 1/9 KI ^ 3 . 32 ) 

f 4a ^ ^ g ±.)^s 

^ ^ -— ) (l-rt(t-t ) ) 

S. s„ o 


’1 


oo . ut u— bf - /- 

- -1_ / i e"“'= [l-e °H(t-t )] Ho Sj <30 

TTO _-ll L 


Jis^ b£ 


. _1_7 i e"“^ [l-e'"*'°H(t-tJ] sin du 


^^2 b' ^ 




- ( 


e i + S. -) (l-«(t-t ) ) 

Q _ O 


-i- / i e'"’^ [l-e 
"=1 b ' “ 


i e-""^ [l-e“''° H(t-t„)] sln(^)^''^s^ du 


ns 


oo 

/ 


2 b' 


i 

u 


j-ut [i-e H(t-t^)] S 2 du} 

(3«33) 


where = x^+y^+(z“hg) ^ , S 2 = x +y +(z+hg) 

It is noted that equations (3*30) and (3*3i) can be 

obtained by the following relations : 
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::2(x/y#z,t) = / C^(x,y,z,t') df' 


^2 (x#y/Z,t) = o Si (x,y#Zjft') dt' 


(3 .34) 

(3.35) 


and = W_» 
o c 

It is foxind that equations (3*32) and (3»33) can 
also be derived by the following relations : 

t '*-’*'*'0 

C 3 (x^y*z,t) = / Cj_(x>y,z,t')dt'-H(t-t^) J (x,y/ z,t'-tQ)dt' 

o o 

(3.36) 

t t“t 

Cp 2 (x/y/Z*t) =/ Cp^(x,y* z,t' )dt'-H (t-t^) / (y,y,z,t'-t^ 

and Wq = W^. (3^37) 

When kj. = 0.0# = 0»0, kj^ = 0.0/ the expressions for C and 

Cp in each case are same as obtained in Chapter II* 

3.4 RESULTS and DISCUSSION 

To see the effect of various parameters on the 
concentration distribution of both the species in each case, 
we use following dimensionless quantities (As in Chapter II)* 


t = t / X = 

hi “s 


/ y = , Z = JJ- / 


h^ Tc . k h^ _ 

]c = -W * '^g ~ D ' ^p “ \ 


k h^ 
_e s 


D 




D 

U h 


1 1,2 

k h 
a s 

D 


^ k h^ 

-b = --D^ ' ^r = ^ U = -D 

_ h D __ W D 

5 _ S — c ^ ^ - _o — ^ 

p ^c P w 1? 

c s 


h D 

S. c — ^ c 


)dt' 


t 


(3,38) 
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The concentration distribution for each case in 
dimensionless form can be written as (dropping bars for 
convenience) : 


(i) Instantaneous flux 


exp b -h't 


Ci(x,y,z,t) = ”—— 372 "— ^2! ® ^ 


+ e 


"2 

4t- 


3 -b't r 


4t 


+ e 


12 

4t 




W exp(2S) -b't - 

Cp(x^Y,2^t) = ^/2 ^ ^ 


2 

h 

4t 


+ e 


2 

_2 

4t' 


-b't r 
- e ^ [e 


2 
4t 


“2 

4t 


t e 

where m' and n' in dimensionless forms are 


( 3 . 39 ) 


(3.40) 


m' = 


n = 




^ + k + kg * k^ 




- * + h> + kp 


2 2 2 , ,2 

s^ = X +y +(z-l) , s. 


2 2 / ^2 
X ty +(z+l) * 


(ii) Constant flux 

t 

C 2 (x/y, z,t) = J Cj^(x,y,Zy t' ) dt' (3 •41) 

o 

t 

%2 (x,y/Z,t) = i (x,yfZ,t') dt' 


(3*42) 



49 


(iii) Step function type flux 

/ t-t 

C 3 (x^y, 2 ,t) =/ C-(x/y,z,tMdt'-H(t-t ) J ° C- (x,y,z, t'-t )dt' 
o o'*' ° 

(3.43) 

t t-t 

Cp3(x,y,z^t) = J Cp^(x/y,z,tOdt'-H(t-t^) J ° (x,y,2,t'-t^)dt' 

(3.44) 

The dimensionless concentration distributions of C and 
Cp at the ground level have been computed and depicted in figs« 

(3. 1-3. 9) for each case and for different values of 
U = 5 .o< =5 1,0 and k^^ = k^^ = k^ ~ ~ 0.t0* It is found that 

C is greater than Cp when k = 1*0, k^ = 10*0 in each case 
(see figs* 3*1 - 3*3) but reverse is the situation when k^ = 1*0, 
k = 10*0 (see figs* 3*4 - 3*6). However, when k^=k=1.0 (see fig* 

3. 7-3. 9), the concentrations G and C tend to their equal equilibr- 

hr 

ium values as downwind distance increases after particular time* 

To see the effects of removal rate parameters on C and 

Cp/ the expressions given by equations (3*39) - (3*44) are 

computed and plotted in figs. (3*10 - 3.12), for following set 

of parameters ; U = 5*0, W = l.o, k = 1 .0, k^ = 1.0, k^ = 0.005, 

kg = 0*05, kj^ = 0.002 , kp = 0.05, z = 0.0, y = 0*0. It is 

seen from these figs, that the concentrations C and C decrease 

P 

as removal parameters increase. Also from these figs* (3*10 - 
3.12), the following remarks can be made as pointed out in 
Chapter II* 
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(i) When the flux is instantaneous/ it is noted from 
fig* 3#10 that as -cime increases the concentrations 
of both the species decrease and the point of 
maxima in the concentration distance profile of both 
the species moves away from the source* 

(ii) In the case of constant flux (see fig* 3*11) it is 
seen that the concentrations of both the species 
increase as time increases and reach to their 
respective steady state values* 

(iii) For the step function type flux, it is observed that 
the concentrations of both the species decrease as 
time increases for t > t^« However for t < t^ the 
profiles for C and Cp will be similar to the case of 
constant flux# 
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FIG3.3 FLUX IS STEP FUNCTION TYPE AT THE SOURCE 



FIG 3. 4 FLUX IS INSTANTANEOUS AT THE SOURCE 
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FIGS. 9 FLUX IS STEP FUNCTION TYPE AT THE SOURCE 



FIG 3.10 FLUX IS INSTANTANEOUS AT THE SOURCE 
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CHAPTER IV 


DISPERSION OF A REACTIVE POLLUTANT FROM A TIME DEPENDENT 
POINT SOURCE FORMING SECONDARY POLLUTANT : WITHOUT 

INVERSION CONDITION 

4 «1 INTRODUCTION 

As pointed out earlier the dispersion of air pollutant 
in atmosphere from a point source is governed by the process 
of molecular diffusion and convection and depends upon factors 
such as stack height, wind velocity, temperature inversion, dry 
deposition, rainout/wa shout, etc^ (Pasquill, 1962; Seinfeld, 
1975; Dobbins, 1979; Novotny and Chesters, 1991), Various 
investigations have been carried out to understand the 
dispersion of air pollutant by considering some of the above 
mentioned factors (Smith, 1957; Pasquill, 1962; Heines and 
Peters, i973a, 1973b; Lamb and Seinfeld, 1973; Exroak, 1977; 
Peterson and Seinfeld, 1977; Karamchandani and Peters, 1983). 
Effects of removal by chemical reaction, rainout/wa shout etc. 
on the dispersal process have been studied by considering them 
as first order processes (Slinn, 1974, 1930; Scriven and Fisher, 
1975; Nordlund,i975; Sander and Seinfeld, 1976; Prahm et al#, 
1976; McMohan et al,, 1976; Sheih, 1977; Calvert et al,, 1978; 
Peterson and Seinfeld, 1977* Alam and Seinfeld, 1931; Hales, 
1982), In particular, Alam and Seinfeld (1931) have studied 
the dispersion of sulfur dioxide and sulfate from a point source 
by solving the steady state three dimensional diffusion equation 
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the secondary pollutant disperses in the same manner into 
the atmosphere as the primary species^ The wind speed is 
assumed to be sufficiently large so that diffusive transport 
in the wind direction can be neglected in comparison to 
a dv action* 


The unsteady state diffusion equation governing the 
concentration C for the gaseous pollutant can be written as 
follows : 


ac 

3t 


+ U 


ac 

ax - 


9 2 - K 

az - 


a^c 


+ K. 


.2^ 

a c 

az 


- (ktk )c 


(4.1) 


where x axis is taken along the wind direction, a towards 

heiaht, U is the mean wind velocity and K , K are diffusivity 

coefficients in y-, z-directions respectively . The constant k 

is the rate of conversion of primary species to secondary 

pollutant, k is its removal rate by scxne other mechanism 
y 

(say washout) and w is the settling velocity. 

The initial and boundary conditions for (4»l) are 


C(x,y,z,t) =0 at t = 0 x,y,z > 0 (4.2) 

C(x,y,z,t} = 6(y) 6(z-hg) at x = 0 (4.3) 

C(x,y,z,t) =0 as y ■* jH» , t > 0 (4.4) 

K t wc = v^ C at z = 0 , t > O (4.5) 

Z d 2 Cl “ 

C =! 0 as z-*«’t>0 (4.6) 
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where is deposition velocity of the gaseous pollutant on 
the groxindw The boundary condition (4*3) implies that the 
concentration has been prescribed at the source in terms of 
time dependent flux W(t)» The following forms of W(t) are 
considered in the analysis : 

(i) Instantaneous source, 

W(t) = Wq 6(t) 

(ii) Constant source, 

W(t) = (constant) (4*7) 

(iii) St^ fimction type source, 

W(t) = W_ 0 < t < t^ 

= 0 t > t 

o 

where t^ is the duration of release of the gaseous pollutant 
from the source# 


Similarly the differential equation governing the 
concentration Cp of the secondary pollutant can be written 
as 


3C 3C 

at ax 


ac 


K — 2^ + ^ + T<C 

^ ay 9z 


Tc 

P P 


(4.3) 


where k is its rate of removal# 

If there is no direct emission of secondary pollute nt 
from the source then the initial and boundary conditions for 
Cp are 

Cp(x,y, 2 ,t) = 0 at 


t = 0 


(4-9) 



56 


Cp(x/y/ 2 ,t) = 0 
Cp(x^y#2,t) = 0 


s = 0 


at 

X * 0 

t 

> 0 

(4 *10) 

as 

y -f-oo 

t 

> 0 

(4*11) 

'P 

as z = 

0 

t > 0 

(4 *12 ) 






as 2 -*• 

C30 

t > 0 

(4.13) 


where Is deposition velocity of secondaiy pollutant on the 

ground* 

Using the following dimensionless variables 

h_ 




h 


f z 


Uh. 


a 


2 V, 

# ri 

s 

o 


(4*14) 


C = 


U h! 

c 

IT" 


c / a 


u hfc 
— ±J 2 . , 

w ' 

c 


w 


2 

'' \ - Wo *’a .-if. , w(t) 

-r— ' Wq = rTT ' = w 

2 C 2 


K 


equations (4*1) -(4 *13) can be written in the dimensionless form 


as (dropping bars 

for convenience) 



9c . ac 

mmu tutrn *Bi. miwi 4M* 

dt ax 

tr 

W ^ 

= 0 

. a^c 

TT 

ay 

32 

(a + ag)c 

(4 

.15) 

c(x,y, Zft) 

= 0 


at 

t = 0 


(4 

*16 ) 

C(x,y,s,t) 

= w(t) 

6 (y) 

6(s-hjg) 

at X = 0 

(4 

.17) 

C( x,y/ 2/ 1) 

= 0 

as y 

■foo 


(4 

.13) 

-I- ttC - 

a 2 ^ = 

NC 

at 

2 = 0 



(4 

.19) 

0 

ii 

0 


as 

z ■* 00 



(4 

. 20 ) 
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(i) W(t) = 6(t) 


(ii) W(t) = 1 


(4.21) 


(iii) W(t) * 1 


0 < t < t 


t > t 


3c 3c 3c a^c a^c 

at^ + - w = /3 + 

ay az 


a c 
P P 


( 4 * 22 ) 


Cp(x,y,Z/.t) =0 at t = 0 


(4-23) 


Gp(x/y,z,t) =0 at X = 0 


(4.24) 


c^(x,y/Z/t) =0 as y +00 


(4.25) 


4- ^ C = N C at z = O 
az p p p 


(4.26) 


C = 0 as 2 •* 0 ° 

ir 


(4.27) 


where 

K k h^ k h 

2 z ^ z z 


k 


d a 


and h is soma characteristic height chosen as source height- 


4.3 METHOD OP SOLUTION 


Equations (4.15) and (4*22) can be written in the 


compact form as 
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P a+a o 

1 -a" a„ 


= 0 


U.23) 


39 3 3 ^ 3 ^ 

where L = - + w -h ^ — j- + — »- is an operator* 

ay 3z 

Using the same transform (see Astarita et al*, (1975)^ 
Aiam and Seinfeld (1981), Chapter II), we get uncoijpled system 


LC - (a + a ) c = 0 

y 


LB - Op B 


= 0 


aC 


where b = ^ ^ ^ 


p g 

The corresponding conditions for B are. 


(4 -29) 
(4*30) 


B(x,y, 2 ,t) =0 at t = 0 


(4-31) 


B(x,y,z,t) = 


-a 


a a - a 
p g 


W(t) 6(y) d( 2 -hg) at x 


= 0 


B(x,y, 2 ,t) = 0 


as y -I* 00 


(N - N)ac 

t WB = N B + — ^ 

3 ^ P (a-a-a) 


at z = 0 


P 


B = 0 as z ■* <x> 


(4 -32 ) 
(4 .33) 

(4*34) 

(4.35) 


Assuming that the solutions of equations (4.29) and (4 *30) are 
of the following forms 


- - X“ X 


C = e 


(4*36) 


w 


B = B 


- " 
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we have/ 



+ — ■L 

ax 


2 

^2 

it"" 

= 3 

a 

~"r“ 

ay 

a a 

« Jb 

32 




2 

-^2 

9B. 


= 3 

a B^ 

. ^ ®1 

at 

t 

ax 

rT“ 

ay 

[|(N 

1 N 

1 

I 

■f 


(a + ttg) 


a 




(4 ,37) 

( 4 ,38 ) 




(a *39 ) 


= 0 * t 0 y-,+.oo^ZH-oo 

y( z-h ) 

= W(t) 6(y) 6(z-hg) e ^ at X = 0 


(4,40) 


(4,41) 


B 


1 


Ta~ -a ' -S ~T eCz-hg) 

P 9 


e 


w 

I 


(z-hg 


) 

at X 


0 (4*42) 



9z“ 


(N - j)Cj_ at 2 = 0 




C£(W^-N)C. 

■*' Ta-ferr ^ 

p g 


0 » 


(4.43) 

(4.44) 


The solutions for both and are obtained by 
using Laplace and Fourier transforms in each case and finally 


we have 


(i) For instantaneous source 
C^(x,y,Z/t), = Wq P(x/y/z) 6 (t-x) 
B^(x/y/Z,t) = Wq Q(x,y/z) 6 (t-x) 


where 
P (x/y, z) 


exp(~ y'*'/4^x) -(a+ag)x 


-(z-hg)- 


471 X 




4x 


(4,45) 


(4 .46 ) 
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( z+h„)' 


~ (47rx) (n - ^) e 


(N-y)^x+CN- ^)(2+h ) 
2) e 2 2 s 


Zirri •« 

X erfc { — ~ + (N - . 

2 V^x 2 J 




Q(x/YrfZ) = X^nTZ??: 7 “ — 


^ « P 
— e ^ 


^ X -(a+a -a )x^ . rr-r 

[(n-n ) ; e ^ P ^ — {2e 

P o VliP" 

2 2Vx' ^ 


-z 

i ® 4 (x 


(N -^)z + (N -^)^(x-xM 

.(N - 2 ) e P ^ P ^ X 

P 2 


-(a-h_)^ 


erfc (-- -i- (N - f ) (x-x' )^^^)}dx' - { (e 

2/x^ P ^ >/4jix 


-(z+hg)' 


1/9 W (n •-J)^X+ (N -j)(z+h ) 

»i/2 p 2 p 2 S 


) - (4jtx)'^'^ (N - j)< 


erfc (• 


2 Vx 




and Dirac delta function 6 (t-x) is defined as (see Carslaw 
and Jaeger, 1941) 

6(t-x) =0 t < X 


X £ t < xte 
t > x+e 
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(ii) For constant source 

= P(x/y,z) H(t-x) 

= Q(x/y#z) H(t-x) 

where H(t--x) is Heaviside function defined as (Carslaw 
Jaeger, 1941) 


H(t-x) = 0 

“ xte 
= 1 


t < X 

X < t < Xte 
t > x+e 


(iii) For step function type source 


= P(x/y,z) [ H(t-x)-H(t“t^-x) H(t-t^) ] 

= Q(x,y,z) [H(t~x)“H(t-tQ-x) HCt-t^)] 

It is foimd that as t_ -* «> equations (4»50) and 
reduce to (^*47) and (4.48) respectively-. 

Finally using (4.36), the solutions for C(x,y, 
B(x,y,z,t) in each case can be written as follows : 


(i) For instantaneous source 
C(x,y,z,t) = Wq P(x,y,z) 6(t-x)e 




^(z-hg) - 


w 


(4.47) 

(4.43) 

and 

(4 .49) 

(4 .50) 
(4.51) 
(4.51) 

;,t) and 

(4.51) 


B(x,y,z,t) = Wq Q(x,y,z) 6 (t-x)e 


(4.53) 
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(ii) For constant source 


C(x,y,z,t) = P(x,y,z) H(t-x) e 


W/ , s 

- 2-(2“hg)- 


U *54) 


B(x,y,z^t) = Q(x,y, z)H(t-x)e 


_ W/ * \ w 


X 


(4*55) 


It is noted here that if a = 0*0# a = 0*0 and t •* «> 

g 

equation (4»54) reduces to same form as obtained by Eirmak(l977) « 
(iii) For step function type source 


C(x/y, z, t) 


= P(x/y, z) [H{t-x)-H(t-tQ-x) 



)- 



HCt-t^) ] 


(4 *56 ) 


B(x,y,z^t) = Q(x#y,z) [H(t-'x)-H(t-tQ**x)H(t-t^) 




(4.57) 


The dimensionless concentration C can be obtained 

Jl 

from C and B in each of the above mentioned cases by the 
following relation 

ac(x/ Yf z, t) 

C^(5c,y/Z,t) = B(x#y/Zt) t (4*58) 

(ap-a-ag) 

4.4 CROSS WIND INTEGRATED SOLUTIONS FOR LINE SOURCE 

In the case of line source # C and C can be obtained 

Ir 

from the above solution in each case separately by integrating 
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C and C with respect to y between -o® to «»• These are given 

XT 

as follows s 

(i) For instantaneous source 

C(x/Z,t) = P(x,z) 6 (t-x) ( 4 » 59 ) 

ac(x, z, t) 

C (x^z,t) = B(x,z,t) + ( 4 * 60 ) 

p g 

where 


B(x/Zyt) = Q(x^z) 6(t-x) 


2 -(z-h -(z+h 

-(a+a +. Si--)x- -Cz-h ) ^ ® — 

g 4 '" 2^"'=' [(e 4X 


+ Q 4x j 


'^47IX 


(n~^ )^x + (N“ y) (z+h_) z+h_ _ 1/0 

•( 47 Ix)^^^(n- “) e ^ ^ erfcf (N- ^)x ^ }] 

2 2 Vx ^ 


Q(x,z) = X^r^:^:=ar7 ■** “ f^^-hg)} 

' p g ^ 

2 


-h 


^ X “(a-HX -a )x' , err? 

[(N-N ) fa 9 P -I— f2 

P o '/^47TX^ 


•( 47 Ix' (N - ^) exp{(N — x' + (N — y)h } 


2 s 
2 


h 


-z 


X erfc( — '^ x)} { 


W\ xrrzT\ •> «■ ■ 1 ® 4 ( X“'X‘^ ) 


2 V x’’ '/7 t(x“x'’ 

(N_ - Z+(N - y)^ (x-x^ ) 


, 'J - 55-; Z+UN - y 


erfc(- 


2^" x-x*^ 


+■ (N„ - y) Vx-x' ) } dx' 
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-(Z-hg)^ -Cz+hg)^ 

I f® 43c ® 4x' 


(47Ix)'''^^ (N - ?) X 
P 2 


exp{(Np-|)2x + (Np~ |3 (z+hg)} erfc( ~— + (N - H) /x) JJ 


( z+h_) 


2 /x 


P 


(ii) For constant source 

C(x,Z/t) = P(x,z) H(t-x) 

, . oc C(x#z#t) 

C (x,Z/t) = QCxvz) H(t-'x) + 

(Op-a-ap) 

(ill) For step function type source 
C(x,Z/t) = P(x/z) [H(t-x) - H(t-tQ-x) H(t«-t^)] 

Cp(x/Z,t) = Q(x/z) [H(t-x) - H(t-t^-x) HCt-t^)] 

a C(x,Z/t) 

+• * 

(cx —a a ) 

'• p g 

4.5 RESULTS AND DISCUSSION 


(4.61) 

(4.62) 


(4.63) 


(4*64) 


The effect of various parameters on the 
steady state concentrations of both species along the central 
line (0/O,l) are shown in figs<* (4<*1 - 4*9) in all the cases 
and for different valuesof t^ W = 1*0. N = 0.4# N = 0.04# 

ir 

z = 1.0# hg = 1.0# 0=1.0 and y = 0.0. It is observed that 
the concentrations C and decrease as downwind distance 
increases (see fig. 4.3). As a increases the concentration 



65 


C decreases but increases in all cases at a particular time 
and location (see fig* no# 4*7)* From these figures it is 
also noted that the effect of r^oval is to decrease the 
concentrations of both the species# 

When the source strength is instantaneous the 
concentration distributions along the central line are shown 
in figs* (4*1 -■ 4 #2 ) for different valuesof t, W^ = l#0, 
a = 0#il^ ttg = 0#02, 0#08 and = 0*2, 0#3» It is noted that 
as time increases both C, C decrease, C being less than 

ir It 

C for all the set of parameters chosen in graphs (4 #1 ~ 4 *2 ) # 

When the source strength is constant, the central 
line concentrations of both species are plotted in figs# (4#3~4#4) 
for different valuesof t, a = 0*11, (Xg. = 0#02, 0#08 and = 0*2, 
0*8, w = 0#04# From these figures, it is observed that the 
concentrations of both the species C, C increase as time 

ir 

increases and reach to their respective steady state values 
as t -• 00 . It is also noted that for t < x the concentrations 

of both the species are zero implying that the pollutant front 

/ 

has not reached the point x as yet* 

When the £Lux is given by step type function, the 
central line concentrations C and decrease as time increases 
for t > t^^ (see figs# 4*5 - 4*6) and for t < t^ the behaviour 
is similar as in the case of constant flux* 

The effect of settling velocity on the unsteady state 
ground level concentration is also studied and it is found that 
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tha concentrations of both primary and secondary species 
increase as settling velocity increases (see fig* 4*9 for 
the case of constant flux)» 

Similar results/ as discussed above/ have also been 
noted in the case of line source for all the three types of 
time dependent sources (graphs are not shown)* 

To see the effect of settling velocity and other 
removal mechanisms (chemical reaction, dry deposition) on the 
concentration distribution at different heights, the steady 
state vertical distributions of the primary species in the case 
of constant flux for both point and line sources are depicted 
in figs* (4-10 “ 4*15)* It is noted fron figs* (4*10, 4*13) 
corresponding to point and line sources, that as settling 
velocity increases, the concentration at a point increases 
when it is located below the source height for a given x* 
However, the results are just the reverse for points located 
above the source height* 

The effects of removal mechanisms on the vertical 
concentration distributions in the case of point and line 
sourcesare depicted in figs* (4*11, 4*12,4*14,4*15) # It is 
noted that the concentration decreases as a increases in both 
the cases* 

The effect of ground level deposition parameter is 
shown in figs* (4 *12,4 *15)* It is observed that as the ground 
level deposition parameter increases the vertical concentration 
decreases, the effect of deposition being more pronounced below 


the source height 
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CHAPTER - V 


DISPERSION OF A REACTIVE AIR POLLUTANT PROM A TIME 
DEPENDENT POINT SOURCE FORMING SECONDARY POLLUTANT: 

WITH INVERSION CONDITION 

5*1 I NTRODUCTION 

In the previous chapter, effects of cheitiical reaction, 
settling velocity, dry and wet depositionson the unsteady state 
dispersion of a reactive air pollutant frcm a time dependent 
source have been discussed (without inversion condition)- In 
this chapter the same problem is studied under the inversion 
condition- In this case#thera exists a layer above the ground 
at height z = H (say) which inhibits the vertical mixing of 
air pollutanti* Mathematically this aspect can be represented 
by 

J^(x,y,z)= g^-wc = Oatz=sH 

where is flux in the z-direction, is diffusion coefficient 
in z-direction, w is settling velocity of the pollutant- 

As in the previous chapter, the following forms of 
flux' are prescribed at the source to study the concentration 
distributions of both primary and secondary pollutants : 


(i) 

instantaneous 

source^ 

(ii) 

constant source. 

(iii ) 

si tsp function 

type source- 
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The case (iii) has also bean applied to study the 
leakage of MIC gas in Bhopal, India* 

5*2 MATHEMATICAL FORMULATION 

Consider the dispersion of a reactive air pollutant 
in presence of an inversion layer fron a time dependent point 
source located at height h^ above the groiond and forming a 
secondary pollutant* The unsteady state diffusion equation 
governing the concentration C of primary pollutant can be 
written as follows, 

where the x—axis is along the wind direction, z—axis towards 
vertical height, U is wind velocity and K^fK^ are diffusivity 
coefficients in y-, z-directions respectively- The constant 
k is the rate of conversion of C to secondary pollutant and k^ 
is its ronoval rate* 

The initial and boundary conditions for C are 
C(x,y, 2 ,t) =0 at t = 0, x,y,z > 0 
C(x,y,z,t) = 6(y) 6(z-hg) at x = 0 

C(x,y,z,t) = 0 asy-»±« t>0 

K^. 'll + wC = v^ C at z = O t > 0 


(5.5) 
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Assuming that there exists an inversion layer at 


height H which inhibits vertical mixing, we have 


+ wC = 0 at z = H 


(5*6) 


As in the previous chapter, the following forms of W(t) 
are considered in the subsequent analysis : 

(i) Instantaneous source 


(ii) 


w(t) = 6 (t) 

Constant source 


(5*7) 


(iii ) 


W(t) = (constant) 

Step function type source 


W(t) = W, 


0 < t < t^ 


t > to* 


The atmospheric diffusion equation governing the 
concentration C of secondary pollutant can similarly be written 

It 


ac ac ac 

3t ^ ^ 9x az 


2 2 ^ 

d^c a c 

K — ^ t K — t kC - k C 

y ay az 


(5*8) 


where 3<p is the rate of removal of G^. 


If there is no direct emission of secondary pollutant 
from the source the initial and botindary conditions for Cp can 


be considered as 

Cp(x,y,z,t) = 0 
Cp(x,y,z,t) = 0 


at t = 0 


at X = O 


(5*9) 


(5*10) 
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C (x,y/Z/t) =0 ss y •* + 

Jr* “*** 


ac 

K ——2. 4" w 
Z 0Z 

ac 

K„ + W =S 0 

Z 3Z p 


Cp = Cp at z = 0 


at z = H 


(5.11) 

(5.12) 

(5.13) 


Using the following dimensionless variables 


K t K X 

r z **** z **** 
t = 1 X = — y / y 


y 

H 


2 = 2 
/ z - H / 


H 


C = 


UH^C 


C = 


UH^C 

I 

w 


/ w = 


K 


W(t) = , Wp 


W 2 

-£ Si 

W K 
c z 


(5.14) 


equations (5.1 - 5*13) can be transformed in the dimensionless 
form as (dropping bars for convenience) 


iS + ^ 

3t 0X 

C(x/y/ z, t) 

C(x/y/Z/t) 
C( x,Yj Zjt) 


ac p a^c , a^c 
ay az 


(a+a„)c 

g 


= 0 at t = 0 


W(t) 6(y) 6(z-h^) at X = 0 

s 


= 0 as y -► + 


t wC = NC at z = 0 

0Z 

■|~ + wC = 0 atz = l 
d ^ 

(i) w(t) = 6(t) 

(ii) W(t) = 1 

(iii) W(t) =1 0 < ^ < % 

t > t^ 


(5.15) 

(5.16) 

(5.19) 
(5.18) 
(5*19) 

(5.20) 

(5*21) 


0 
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3 C 3C 3C^ 

3^C 

— 3 — E. ^ 

3Y^ 

9^C 

^ + ac - 0 ^ Cp 

(5*22) 

Cp(x,y,z,t) = 0 

at t = 0 


(5.23) 

Gp(x,y^z,t) = 0 

at X = 0 


(5.24 ) 

Cp(x,y,z,t) = 0 

as y -» + 

oo 

(5.25) 

3C 

„ Cp = Np Cp 

at z = 0 


(5.26) 

+ w Cp = 0 

at z = 1 


(5»27) 

K 

where ^ ^ , a =■ 

^z 

kH^ 

k v^H 

— M ^ 



K ^ ^ “ K 

z z 



V' 

QL SS aiMinkMaMMiw 


V. 


H 


/ Np 


K 


. 


z ^ z 

5.3 METHOD OF SOLUTION 


As in Chapter II, following Astarita et al. Cl979), 
Alam and Seinfeld (1981), equations (5«15) and (5*22) can be 
written as uncoupled system j 


LC - (a + ag)G = 0 

LB - a B = 0 

p 


where 
L = - 


B = Cp(x,y,z,t) 


3y 

ac(x,y,z,t) 
(Op - a-ttg) 


3 -3 ,3 n 3“ , 3^ 


(5«28) 

(5.29) 


(5*30) 
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The boiindary conditions for B are 
B(xyy/Z,t) =0 at t = 0 


BCx/Y/Z^t) = 


-a 


a -a-a 
P g 


W(t) 6 (y) 6(z-h^) at X = 0 

S 


BixtYfZft) =0 as Y •* ± 


a(N -n)c 

= (Np-w)B +• x’ap-cx-dgT 2 = 0 


+ wB = 0 at z = 1* 
a ^ 


Introducing the transformations 
2 


w 


w 


C = e 




B = Bj e 




we have from equations (5 #2 3 - 5»35) 


at 

ac- 

+ i 

ax 

= e 

ay 

■ t — •“ 

aB. 

as^ 

2 

3 


3t“ 


_ |3 — j. 

By 

+ - ( 


P 1 


and boundary conditions are 


C^ = 0 asY’*i°°# att = 0 

B^ = 0 


(5*31) 

(5.32) 

(5*33) 

(5*34) 

(5.35) 

(5*36 ) 


(5.37) 

(5.38) 

(5.39) 
(5*40) 
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w 


= W(t) 6(y) 6(z-h^) 


at X = 0 


■^Ct y(z~h ) 

*<y) 6 ( 2 -h ) ® 

p p 


at X » 0 


!£? 

az' 


= ~ 2’^^! z = 0 


ac 

3 z 


^ "** z ^ z = 1 


9®1 w a(N -n) 

~ = - j)B^ + 5 j~=2_ c ^ at z = 0 

p g 


aB- 

aV * 2 ®1 = 0 at 3 = 1. 


(5*41) 

(5-42 ) 

(5.43) 

(5.44) 

(5-45) 

(5.46) 


The solutions for both and are obtained by using 
Laplace/ Fourier transforms and method of separation of 
variables* We have 


= Wq P^(x/y/z) 6 (t-x) 
~ ^o 6 (t-x) 


'-A’-here, 


exp (-y^/43x) 
P’"(x/y/z) ZZri 


2 




X s 
n=l 


cos sin y) 


© 


(5.47) 

(5.48) 
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expC-y'^APx) ^ 

Q''(x/y/z) = ' (a -a-oc } 

'^47rFx p g 


oo 

X S 


2 <l^n + l-> (''n sin 


“(a + Ji^)x 


(N -N ) (2+w(i~z)) 


X { 


p .. . . 

( w + N ^ ( 2 +w ) ) 
pi 

— (ji^ sin Mn + 


2 


(?v. COS X^h -fN- 
X ^ ms 1 

sin 


m=l 


(X^4-N^) (A^ + ^ + f^+Nj_(?^ + ) 


/ e 
o 


“g p Tn m 


4 


- cos Mn hg + sin Mj, h^)} 


w 


N- = N - o 

X ‘O 


w 


N , = M - ^ 

pi P 2 

and 7 , 's, u 's are roots of following transcedental equations^ 

X (N. + y) 
n 1 2 • 

tan X = — 5 “ 

Wi 


‘‘n = TTJ w 




and 6(t-x) is Dirac-delta function, 
(ii) For constant source 


= P*(x,y(rz) H(t-x) 


dx' 


(5.49) 
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= Q*(x,y/z) H(t-x). 


(5*50) 


It is observed from above equations that and 
are independent of time as t -♦ «>♦ This implies that both 
concentrations C and C reach to steady state* 

ir 


(iii) For step function type source 

= p'^(x,y/z) [H(t“x) -H(t-t^-x) H(t--t^) 1 

= Q*(x,y/z) [H(t-x)*-H(t-t^-x) H(t~t^)] 


Finally the solutions for C and B in each case 
can be obtained as follows : 

(i) For instantaneous source 


C(x,y,z,t)=WQP’‘‘(x^y/z)6(t-x)exp{- j(z~hg) - ~x} 

^2 

B(x.y^z,t)=WQQ*(x,yrz)6(t-x)ej^£<- 


(ii) For the constant source 

2 

C(x.y.z,t) = P*(x.y/2) H(t-x)exp {- j (z-h^) - —Jci (5.55) 

B(x,y,z.t) = Q^^x.y,z) H(t-x)e:<p{- | (z-h^) - * <5.56) 

It is noted here that as t ■* «> and w = 0 the above 
equations reduce to the same form as obtained by Alam 
Sein f eld (19 81 ) . 


(iii.) For stop function typo source 
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C(x/y,z,t) = P*(x,y,z) [H(t-x)-H(t-t -x) H(t-t )] x 

o o 

X exp £~ j(z-hg) - ~x} (5.57) 

B(x/y/z^t) = Q*(x,y,z) [H(t~x)-H(t-t ~x) H(t-t^)l X 

Q O 

2 

X exp {- ^(z-hg) - -^-x}* (5»53) 

5.4 SOLUTIONS FOR THE LINE SOURCE 

The solutions are obtained for the line source by 
integrating equations (5.53 - 5.58) with respect to y between 
-oo to oo and are written in each case as follows ; 


(i) For instantaneous source 


C(x|Z/t) =s P^(x. z) 6(t-x) 

C (x,z,t) = W Q*(x,z) 6 (t-x) + 

where 

J2 

P^(x/z) = exp {-(ct+ag + ~ ^(z-hg)} 

2 

«> 2 (X.^ ^) (X cos X h + sin X h ) 

X S n 4 n ns i ^ 




(5,59) 


(5-60) 


(Aj^ cos Xj^z + sin Xjjz) e 

2 

,, 00 cos Mn^+N 1 sin u z) 

P “gn=i 


2^ , (N-N^) (2+w(l-z) ) 

exp( -(a )x){ 

P " (wfNpj^ (2+w)) 


(4^ sin 4j^+Npj^(i-cos 4jj)) 
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X r 

m=l 


2 

2 hg -f sin X^h^) 

+ ^T» 


X 

X J‘ exp 
o 




2 2 

-a„-^x„ + dx'-(Mj, cos jLi^h^+Np^ sin ju^h^)} 


<3 P 


(ii) For constant source 


C(x^z#t) = P^(x/z) H(t~x) (5.61) 

a c(x/Z/t) 

C (x,z,t) = Q‘’^(x.z) H(t-x) + (5.62) 

(ap-ct-dg) 

(iii) For step function type source 


C(x/Z,t) =5 P*(x/z) [H(t*-x)~H(t~t^-x) H(t-t^)] 

a C(x/Z^t) 

C (x,Z,t) = Q‘’^(x,z) [H(t-x)-H(t-t -x) H(t-t )] t 

P ° ^(a~a-a 

p g 


5.5 RESULTS AND DISCUSSION 


The effects of various parameters on the central 
line (0,0,h„) concentrations of both species in all the three 
cases are depicted in figs* (5.1 - 5.7) for different values of 
t, (3 = 10.0^ N = 2.0. Np ~ 0*2. w = 0*2. y = 0.0, h^ = 0.2. 

It is noted that the concentration C decreases as a increases 
but C increases in all cases at a particular time and location 

ir 

(see fig. 5.7 for exart^le in the case of constant source 
strength). It is also seen that the concentrations of both the 
species decrease as downwind distance increases* The effect 
of removal is to decrease the concentrations of both species. 
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When the source strength is instantaneous the central 

line concentration distributions are shown in figs* {5«1 - 5.2) 

for different values of t^ W = 1*0, a = 2*77, a = 0.0, 2.0 

and a = 0.0, 20.0. It is seen that C and C decrease 
P P 

as time increases, the concentration being lass than C 
for the set of parameters chosen in graphs (5.1 - 5.2)* 


When the source strength is constant, the central 

line concentrations of both the species are plotted in figs* 

(5-3 “ 5-4) for different values of t, a = 2*77, = 0*0, 2.0 

and a = 0*0, 20*0, W = 0-2- From these figures, it is found 
P 

that the concentrations of both the species C, Cp are zero 
for t < x- It inplies that the pollutant front has not 
reached the point x as yet* The concentrations C and 
increase as time increases and reach to their respective 
steady state values as t 


When the flux is given by the step-function type, 
the central line concentrations C and decrease as time 
increases for t > t^ (see fig. 5-5 - 5-6) and for t < t^, the 
behaviour is similar as in the case of constant flux* 


The effect of settling velocity on the unsteady state 
ground level concentration is also studied and it is found that 
the concentration of primary species increases as settling 
velocity increases (see fig* 5,8 for the case of constant flux), 

In the case of line source, similar results have 

been found in all the three cases (graphs are not depicted)- 
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The effect of settling velocity# chemical reaction and 
dry deposition on the vertical concentration distributions of 
primary species in the case of constant flux for both point 
and line sources are depicted in figs. (5,9 - 5,14). It is 
seen frcw figs. (5,9# 5,12) that as settling velocity increases 
the concentration C increases below the source height# but 
reverse is the case above the source height. The effect of 
chemical reaction is to decrease the concentration of pollutants 
in the case of both point and line sources (see figs. 5.10# 5.13), 
The effect of ground level deposition parameter is shown in 
figs. (5.11, 5 * 14 ). It is observed that as deposition at the 
ground level increases the concentration decreases and this 
effect is more predominant below the source height. 

5.6 EFFECT OF INVERSION LAYER ON THE DISPERSION OF POLLUTANT 

To see the effect of inversion layer on the dispersion 

of a pollutant, we compare the corresponding concentration 

distributions in the case of a continuous point source as 

obtained in Chapter IV (without inversion condition) and in 

Chapter V (with inversion condition) , The concentration of 

primary pollutant in dimensionless form (the dimensionless 

quantities is defined in Chapter IV, here h^ is taken source 

height h„) is given as follows in the absence of inversion, 

2 

exp(-y^/4|3x) ~ j(z-l)- (-^ + a + 

C(x,y,z) = — — > e 
4Trx V 0 

-(z-l )^ -■( ztl )^ 

X { [e +e ^ ] - (47rx)^ ^ (N - |) X 
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5.7 application OP CASE III • MIC LEAKAGE IN BHOPAL, INDIA 

The dispersion of Methyl isocynate (Mic) leakage from 
Union Carbide Factory in Bhopal took place on Dec. 3, 1934 in 
the middle of night. It has been reported that at the time of 
episode a light breeze (i.2 kiti/hr) was blowing in south*“so\ith 
ea'=-*-'^riy direction which changed its direction towards south 
south-west and slowed down later. The temperature was around 
12— 14°C at that night. The MIC gas (approximately 40 tons) 
escaped through the nozzle of the 33 meter high vantline and 
leaked for about one hour (A report, 1984)* This gas, being 
heavier than air, settled down over considerable areas of the 
city when the wind died do\>m. The limit of region where people 
got seriously affected was about 8-9 Kms from the factory. 

The report also suggested that the gas travelled over 15 
kilometers covering area of over 50 square Kms, affecting nearly 
two lakh people, constituting one fourth of the city population. 

MIC is strongly toxic and highly unstable* It decomposes 
easily in the presence of moisture forming methyl amine 
according to following reaction 

OH 

ngo + H„0 -* CH. 2 N = C -* CH^NH„ + CO2 

(mic) (Methyl amine) 

MIC affects all living beings including vegetation. 

In a report, it has been recorded that when four persons were 
exposed for 1 to 5 minutes, the following types of effect were 


noticed : 



S5 

0 »4 ppm - no effects 

2»0 ppm - lachrymation, irritation in nose and 

throat 

4 *0 ppm “• Symptoms of irritation more marked 

21.0 ppm - Unbearable irritation of eyas/nose 

and throat# 

In cases of prolonged eiqposure MIC acts by destroying 
protiens and lipids in the lungs# This leads to changes in the 
permeability of the lung membranes which get filled up with 
water causing death * The maximxam permissible ejiposure limit 
for industrial workers prescribed is only one-fifth of one 
ppm for MIC* 

The complete set of data is not available for MIC 
dispersion at the time of episode/ the following 
parameters are^ therefore^ chosen in the computational analysis 
of Case (iii) being applied to Bhopal gas leakage (see Table i) 
(A report, 1984). 
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Using data in Table 1, the expression for C given by 
equation (5.57) is calciilated for different values of w = 0.0,4 .0 
by choosing the effective height of the source as 40 meters. As 
pointed out earlier, the leakage time t^ is taken as one hour 
(i.e. dimensionless t^ = O.OIS). 

The concentration distribution at t = 0*018 (one hour) 
for different values of w = 0.0, 4*0 are shown in fig. 5.17. 

It is seen from this graph that 2 ppm concentration line 
(dimensionless concentration C = 1.4537) intersects the 
concentration profile for w = 0.0 at a point which corresponds 
to a dimensionless distance 0.0163. This is equivalent to about 
10.5 Kms. in the downwind direction* Also the 2 ppm concentration 
line intersects the concentration profile for w = 4.0 at the 
dimensionless distance 0*0172, which is equivalent to 11.5 
kilometers. Similarly the concentration distributions at 
t = 0*02 7 (l|- hour) for different w = 0.0, w = 4*0 are depicted 
in fig* 5*18. It is also observed from these graphs (5.i7*-5.18) 
that as settling velocity increases the ground level concentration 
increases. 

Since the limit for actual affected zone along the wind 
direction in the city (wind changed its direction during the 
episode) varied from 8"’9 Kms. at different locations, the present 
analysis though approximate (wind velocity and diffusivities are 
assumed constant, change of wind direction ignored) provides 
useful information regarding the impact assessment of the 
dispersion of MIC leakage in Bhopal, India. 
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CHAPTER VT 


EPFECT OF VAJUABLE WIND VELOCITY AND DIPFUSI- 
VITIES ON DISPERSION OP AIR POLLUTANT FROM 
POINT AND LINE SOURCES 

6.1 INTRODUCTION 

The unsteady dispersion of air pollutant from a time 
dependent point source under going first order chemical 
reaction forming secondary pollutant has been discussed 
in chapters II-V, The analytical solution of unsteady 
state three dimensional diffusion equation has been obtained 
by taking into account first order chemical kinetics^ dry 
and wet depositionsin chapters IV and V. It is noted here 
that in these studies the wind velocity and diffusion 
coefficients have been taken as constant. In general# 
however# the wind velocity and diffusivities are functions 
of space coordinates and time, but in the case when horizontal 
homogenity is assumed, the wind speed and diffusion coeffi- 
cients may he taken as function of height only. In a e 
cases when the wind speed and diffusion coefficients are 
power law functions of height, the aoalytiCca solution of 
steady state two dlmenslcnal diffusion eguation for non- 
reactivo pollutant has been obtained when the ground surfaoe 
is acting as a roaectlng plane (smith, 1957, Demuth, 1978, 
Robson 1983) . For general funotional forms of wind velocity 
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and di f fusivi ti. ss no ©xact solution of thjTGs dirnsnsional 
diffusion ©c 3 .u,ntion with chemical and othen iremoval 
mechanisms is known. Since a function can be 

approximated by a set of step functions, one way to study 
this kind of problems is to divide the atmospheric layer 
into several layers such that the diffusion coefficients 
and wind velocity are taken constants but different in 
different layers. The diffusion equation is then solved 
in each layer by using these constant values with appropriate 
matching conditions at interfaces of these hypothetical 
1 ayers , 

In view of this, in this chapter, the steady state 
dispersion of reactive air pollutant emitted from point 
and line sources has been discussed by taking into account 
dry deposition and under the chemical reaction inversion 
condition and dividing the inversion layer into two, and 
three layers. 

6. 2 BASIC EQUATIONS AND SOLUTIONS 

(i) Tv^o Layered model 
( a) Point Source 

Consider the dispersion of a reactive air pollutant 
from a continuous point source located at height h^ above 
the ground in presence of an inversion layer (0 ^ z _< H) of 
height H divided into two layers: I (0 < z < hg) and 
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II ~ ^ ^ / ^g being height of lower layer. When 

x-axis is taken along the wind direction, the steady state 
diffusion equation governing the concentration of the 


pollutant in 

each region, can be written as follows 


Region 

I (0 < z < 



ac^ 

^1 ax“ 

a^c^ 
= "^zi az 2 

" ay2 ■” ’"''l 

( 6 . 1 ) 

and boundary 

conditions 

are 


(i) 

at X = 0 

* 0 

(6.2) 

(ii) 

C^= 0 

as y >► + 00 

(6,3) 

(iii) 

ac. 

K = 'V 

zl 3 z 

’’^C^ at z = 0 

(6.4) 

(iv) 

ac^ 

ac 

22 ^ ’ ‘= 1 = S ^ “ "g 

(6.5) 


Region II (hg < z < H) 

ac^ 


u 


2 ax 


d^c a^c 

K . — ? + Ko — 9 


kc. 


a 2 


ay 


and boundary conditions are 

(i) C = ^ 6 (y) 6 ( 2 -hg) at x=0 

(ii) C^= 0 as y -* + °° 

3C 

(iii) = 0 at z=H 

ac, ^ 

^z2 ' ^ 2 " 


at z 



( 6 , 6 ) 

(6.7) 

( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 
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where are the concentration of pollutant in the 

two regions^ are the mean velocities# and 

K^j^/ ^22 diffusion coefficients in y- and z-directions 

respectively, is deposition velocity on the ground, 

k is chemical reaction rate and 6(,) is Dirac delta 
fianction and W is source strength. 

If the source lies in the first region, the boundary 
condition at the source for equation (6.1) and (6.6) are 


W 

u- 


Cj - 0 


at 

X = 0 

(6.11) 

at 

X = 0 

(6.12) 


Using following dimensionless quantities 


K X 

zmax r _ z 

X = 1 , 2 = ^ , 

u 

max 


^i 


^1 V - Z 

u ' ^ “ H 
max 


„ Ta _ h 

T/ — "U 'h, = — — 

K ^ Yj - - , n^ . n 


K, 


zi K 


zmax 


'd K 


zmax 


^g H ^ ^s H 


u 

rn3Xi /n *1 22 1 o^ 

Ci = — =1 - 1-2) 


equations (6,1 - 6.12) become 


Region I (0 < z < h ) 

y 


(6. 13) 


ac. 

= Q 4* V 

ax ay^ 


a^c^ 

1 gTF - “i^^i 


(6.14) 


0 at X = 0 


(6.15) 
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= 0 as y -» (6.16) 


ac. 

— = = N.C. at 2 = 0 

az 


(6.17) 

_ ac-i _ a c,^ _ _ 

K — i = K ^ C = C at 

2^35 ^2g- 12 

" =^g 

(6.18) 

Region II (h^ < i < 1) 



ac„ a^c, a^c 

-^ - e, -rr * ’>2—^' h'^a 

ax ay ^ ^ 


(6.19) 

c = ^ 6(y) 6 (Z'-h ) at X = 0 
^ "2 

j 

(6. 20) 

= 0 y *♦ ,.t ®° 


(6.21) 

- — 2 =s 0 at z = 1 


(6.22) 

3 z 



K 2 ^ = ^zl^ ^ ^ ^2 

ZZ 92 32 

Z = hg 

(6, 23) 


where 


^z2 a * — a = ^ 

r~ ' ^2 = r" ' ^ u ' 2 ^ 

u. ^2 1 2 


K 


^ ^yi 0 =s — ^ 

N. = 2- ^ 

K„ zmax 2” zraax 


a 


kH 

K, 


zl 
2 


il. 




zraax 


/ly a ® raox ("a^ / r. — ) . 

K » raox (K # -K-g 2^^ ' max ^ ^ 

zmax 2 ^ 
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The solutions of equations (6.14) and (6.19) are 
obtained by using Fourier transform^^ with the assurrption 
K ^ K 

/3i= i.e. ^ (dropping bars for convenience) 


1 ^2 
2 


as 


N 


e 


CAx,y,z) = 


(y74pTx) ^2 

1 CO -6_X u, , 

S e " R 


1 

sin S-j -] ^^) 


lln 


n=l 


h 

‘37^ ailn*'g'"'=°® ®UnV 


(6. 24) 


C 2 (x, y, 2 ) 


(y V4(3^x) 


'/'4ti3^x 


S (tan a^^sin 

n=l ^ 


(tan ai2„stn=i2nV'=°®^3f*g' 


(6.25) 


where 


R„ 


(tan sin h^t cos 


n 


PCtan ai 3 „sln h^t cos 


U-P 


r 1 


n: 


n: 


4 a 


[2%i„ (1-^ ^ (1- “T ^ 


lln ^ ^ "!ln 


Uln 


2N. 


Uln 


(cos 2 a^^^hg-l)] 




4oi 


i 2n 


(1+ tan ^i2n^ 


+ (l-tan\^) (sin ^a^^^h^) 

- 2 tan 2 a^ 2 n''‘^°® ^^12n ^g^ ^ 
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itl - a, 

^lln i 


/ 6^ - Ct 

ri " ®l2n “”7; 


P -1 - !/(::; — sin 


^lln "g 


®lln V =“ ®Un V 


Pj - l/(tan ai 2 „ aln cos a^j^hg) 


and 6j^s are the solutions of following transcedental 
equation: 


K 


cos a^ 3 _h^-aj,^ sin a^^h^) 

zl 

(~ sin a^^hg+ cos a-, ^hj 


"11 


"11"V 


‘^z2^H9. ^12 V 

(tan a^ 2 ®^^ a 2^2 ^12 


(6.26) 


NOW/ if the source lies in first region/ the concen- 
tration of pollutant in different regions is given by 

.2 N, 




X 


Cjj^Cx/ Y/Z) 


CO -a^x 


( sin aaa^z+cos aii„z) 

“■fin 


— — 2 e “ R ^ 

n=t {— ^^slnajj^h^+oos h^) 

(6.27) 




e 


1 00 


ViF^Ttx n«l 


6^ (tan aj sin aj a+cos a^ a) 
2 e K 


no 


(tan ai 2 „ ai" ^12nV‘'°® ®12n‘'g> 


( 6 . 28 ) 
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where 


R 


no 


®lln ®Un 

P(tan sin cos 


nnd P is dofjned above. 


(b) Infinite Line Source 


Consider the dispersion of a reactive air pollutant 
from a continuous infinite line source in the two layered 
inversion region as before. In this case, the pairtial 
differential equations governing the concentration of 
pollutant in each region are. 


Region I iO < z < h^) 


u 


3C. 

1 = K i - kC^ 

1 o^c ^^az^ ^ 


with boundary conditions 


(i) = 0 at X ~ 0 

3 C 

(ii) K - = vX, at a = 0 

vxxy rv.^2. gz d 1 




ac 


(iii) -ft == ^z2 “3 


, C- - c at z»h 

d * I y. 


g 


(6. 29) 


( 6 , 30 ) 
(6.31) 


(6,3 2) 


Region II (h < 2 < H) 

9 

3 0 3^ 

'"2 “’^^2 3^2 




(6,33) 
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(i) 


(ii) 


W 

~ 6(2-hg) at X = 0 


3C 


32 


2 = 


0 at z = H 


3c ac 

(ill) az ~ ^2l Tz ^2“ at 2 = h 


(6.34) 

(6.35) 

(6.36) 


If the source lies in the first layer, the boiindary 
conditions at the source for equations (6.29) and ( 6 . 33 ) are 

W 


Cf = ~ 5 ( 2 -h^) at X = 0 


C 2 = 0 


at X = 0 . 


(6.37) 

(6.38) 


The solution of equations (6,29) and (6.33) subjected 
to above boundary conditions can be obtained as before and 
written in the dimensionless form, for each region, as follows 


^1 

CO -el ailn^ + ®Un^> 

c^(x, 2 ) = S R^e ^ (6.39) 


n=l 


( 


a 


lln 


sin an„ h + cos a^i-hJ 
lln g lln g 


0 

00 "^<5 „x (tan ^ sin a^ ^2 + cos a^ ^z) 

C (x,2) = S R e ^ ^ 

^ n=l ” ^l2n ^l2n^g'^*^°® ®’l2n ^g^ 

(6.40) 

When the source lies in the first region, the 
corresponding dimensionless distributions ofC^(x, z) and 
C 2 (x, 2 ) are given as follows: 
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Ni 

C^(x,z) = S R^^e -^- 


(6.42) 


C^Cx^z) 


“ jj j5^ (tan a^^cos a^^zi-sln aj^^z) 
n=l ”° (tan aj^jjjSln aj^^hg+oos 


where various parameters are same as (iefined in the case 
of point source. 

It is also noted that equations (6, 39) -(6, 42) can be 
obtained by integrating equations (6, 24)~(6. 28) with respect 
to y between " and respectively. 


(ii) THREE LAYERED MODEL 
(a) Point Source 

consider the dispersion of a reactive air pollutant 
from a continuous point source of strength W, located at 
height h^ above the ground in the second region of the inversion 
layer which is divided into three layers: I (0 < z i hj^) , 

II (h, < z 5 hj) and III Ch, < z < H). The steady state 
diffusion equation governing the concentration of the 
pollutant in each region can be written as 


Region I (0 < z < h^) 


3 C. 

„ i ta Ic. 

ax 


— h ^zi 
yl 3y2 zl 


d 

dz^ 


- k c. 


(6.43) 


with boundary conditions 
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(i) 

0 

It 

r -4 

0 

at X = 0 

( 6 . 44 ) 

(ii) 

= 0 

at y .» + 00 

( 6 , 45 ) 

(iii) 

ac- 

K 

zl 32 

~ ''^d^l ^"^2 = 0 

( 6 , 46 ) 

(iv) 

ac- 

!/■ -^■.- 

zl 3 2 

^22 Tz ' ^1“ ^2 at 2 = h^ 

( 6 , 47 ) 

Region II (h^ < 

z < h^) 



ac 3 2c 3 2c 

^^2 It " ^^2 772 ^z2 T;! " ^"^2 

( 6 . 48 ) 

with 

boundary conditions 


Ci) 

p „ iL 
^ “2 

6(y)6(2-'hg) at X = 0 

( 6 . 49 ) 

(ii) 

= 2= ° 

as y •* + oo 

(6, 50 ) 


ac„ 3C. 

ac, ac 

^Z2 irl = ^^3 -ai ^ S = S ^ “ ^2 (6.52) 


Region III (h^ < z < H) 
9C 


u 


3 dx 


K, 


3^C 3^C 

y3 ay 2 23 3^ 


kc. 


with boundary conditions 

(i) C„» 0 at X = 0 

J 


( 6 . 53 ) 


( 6 . 54 ) 
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(ii) 

S ’ 

= 0 

as Y ** + <» 

(6.55) 

(iii) 


3C^ 
3 z 

3C., 

= ^2 "at ' S' ^2 =‘=>'2 

(6. 56) 

(iv) 

3C 

3 z 

= 0 

at z = H 

(6.57) 


I£ the source lies in the first region, the boundary 
conditions at the source for equations (6 , 43), (6 , 48) and 
(6.53) are 


C. = 6(u)S(z-h„) at X = 0 

1 u^ s 

C = 0 at X *= 0 (6.58) 

C- = 0 at X » 0 

3 

Using following dimensionless quantities 


h - 2 a S u. * -i. , y = ^ 

X — / 2 TJ/ Uj^ # j: 


^ax 


K, 


K 

■gi" K 


zi 
zmax 


' ^d ■ K, 


v,H „ h, ^ h 

d r as -=• h « — 

' ^i H ' H 


zmax 


u„ H‘ 


(6,59) 


C « C, (i=l, 2,3) 

i W i ' ^ ^ 

equations (6. 43) - (6. 57) can be written, in the cilmensionless 

form as 

Region l (0 < z < h]_) 

.2s 




3 X 


— ^2 * \ 

^ dY 


3 z 


^ - 0-1 Ci 

2 1 i 


(6.60) 
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^1 

= 0 

at 

X = 

0 

(6.61) 


= 0 

as 

y •* 

+ 00 

(6.62) 

Ifl , 






a 2 

= 

at 

2 == 

0 

(6.63) 


ac- 


9C^ 



^2l 

1 _ 

32 

22 

2 

ai 

' ^1 ~ ^2 ^ “ ^1 

(6.64) 

Region II 

< 

2 < 

^2^ 



®=2 
a X 

n 

C 

= 0 - 
2 

Owm* 

ay2 

+ r. 

2- 
a c 

^ a c 

1 9^2 22 

(6.65) 


6 (z~hJ at X = 0 ( 6 , 66 ) 

2 

C^ = 0 as y-» + « (6,67) 


ac^ - 

K ^ =5 K ^ , 

02 ^^3 2 


C^** ^2 ^ ^ ^1 


3C - 3C„ _ „ ^ « 

— * / C„ = C_ at 2 » h, 
a 2 3 z 2 3 


(6,68) 


( 6 . 69 ) 


Region III Ch^ :< 2 < 1 ) 


ac 

SSS rj 

3X 3 

3 y^ 

^'3 • 

3^0 

X - 

ai2 3 3 

(6.70) 

C 3 . 0 

at X = 

0 


(6.71) 

C 3 = 0 

as y •* 

4* 00 


(6.72) 



10 2 


ac , 3c ^ 

jr — 2 = K — — » . C- = n 


■^3 Sz 


z 2 


dz 




(6.73) 


aCg 

az 


at 2=1 


(6. 74) 


where 


K 


zi 


ct. 


u« 

X 


a 




a = 


ml , ^ 

zmax K_ ^ K u . 

zl zmax i 


(i=l, 2,3) 


\2' ^zll %ax= ">ax(u^,U2,U3). 


The solutions of equations (6.60), (6.65) and (6,70) 

K 1 1^2 ^y3 

with the assumption ~ ” u^ ' 

written (dropping bars for convenience) as 


C^(x,y,z) 


C2(x,y,z) 


y ^ 

w’ „ 


2'/%^x n=l 


S «nl 


2fW^ t=l 


2 ^1 


- 

e " G^^i.z) 


** 2 

"■ ” vw 


(6.75) 


(6.76) 




-6’2x 


c^(x,y,z) = T — \l ® ®3n^^^ 


2'/WFp< 


(6.77) 


N 


G^^( 2 ) » i— Sin a^^^^z + cos 


lln 
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^I2n2 + b .[^ z ) 

" ^5 ^l 3 n ^Isn^ •*■ ^os 

P3= l/CN^/a'i^ sin + cos 

P^= 1 /(R Sin a^2j^h^+ cos a][2j^h^) 

P = < 2 n^ 2^^4 

Ctan a-j^3^ sin ^]_2n^2"^ 00s ^2.2n^2^ 

K^2a^2®^^ a2 2^2^^^ a^3Sin a^3h2+cos aj^3h2) 

+K^3COs a;3COs a;3h2-.a‘ 3Sin a‘3h2) 

a^gSin aJ^3h2+cos a[3h2) 

+ iC^sSin a^ 2^1 2! a2^3tan 3^^3003 a^ 2 ^ 2 "'^l 3 ^^^ ^ 13 ^ 2 ^ 


Nr ^''1 , 

)sin 2aJij^h^+2aiin(l+ *7^ )hi 
^lln ^lln 

2Ni , - 

„ (cos J 

^lln 


2 

+ Us [Cl-R^{sin 2a^2n^2“®^^ 2aj^2n^l^‘^2a^2n^^‘*'^^ 

'^^12n 

- 2 R ( cos 2a\ 2n^ 2"'^*^® 2n^l^ ^ 


R, 


nl 


F, 


UiP:r 

. _i_i [d- 


*®Un 
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^ 


+ 2a5^3^(l-.h2) (l+tan^ ^{3^) 

-2 tan a{3j^(cos 2a^3^-cos 


/ 6 ""^-a 
» =v ^ 1 


/a 12 . 


y. ' ®l2n 


/ 6 „ I 5 • ^_cfc 

• =>r-iLr~2. a‘ =:/iiLJi3 


and. 6j^*s ane eigen values of following transcendental 
equation 


I -I 

^z2^l2®^^ ^12^1^’T a^i^j^+cos 


+KaiCos aijh^(H^cos aijh^-a^jsin 


*'■ 22 ^ 12 '^°^ '^12*'!^ I aJj^hj+cos aj^hj^) 

®11 


■K^j^sln aJ^hiCNjCos a’^h^- a^iSin a-^h^) 


a^^^Ctan a(3Sin a*j^3h2+cos a^gh^) 

tjC ^sCOs a2 gC a^gtan a^^nos ^ 13 ^ 2 ^ 

^■^2^1 ^12^2^^®'^ a[3Sin aj_3h2+cos a|3h2) 

-'K^3Sin a|2i^2(a^3tan a23COS ^13^2^ 


(6.78) 
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If the source lies in the first region, the corres- 
ponding dimensionless concentration of pollutant in each 
region is 


X 00 




CT(x,y,z) = — S R' e ^ G, (z) 


(€.79) 


C Cx,y,z) = ' - S e "" "" G^(z) 

^ n=l o 2n 


^ 2 

-‘fe , 2 

O 1 “ "“dr, X 

C Cx,y,z) = — — ^ ^ 

^ 2f%F^ n=l o 


(6.80) 


(6.81) 


P I „ '^In 
^lo 


(b) Infinite Line Source 


In this case the concentration of pollutant, when the 
source lies in second region, is obtained by integrating 
equations (6. 75) ~ (6. 77) with respect to y between - <» to <» 
and written as 


00 -6„ X 


C 3 ^(x,z) = 2 e G 3 _j^(z) 


(6.82) 


ww n 

C (x,z) « S R^l e G 2 „(z) 


(6.83) 


1. 1 2 
-6„ X 


C3(x,z) » S 


? R . e ^ G,„(z) 


(6.84) 
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Similarly, when the source lies in first region, the 
concentration distribution in each region is 




00 




C^Cx, 

z) 

= s 

n=l 

^no 

e Gi„(z) 

(6.85) 



00 





z) 

= E 
n=l 

^o 

e 0 ^^z) 

(6.86) 



CO 




c^Cx, 

z) 

= E 
n=l 

^no 

e 

(6.87) 


6.3 RESULTS AND DISCUSSION 

To study the effect of various parameters on the 
concentration distribution of air pollutant, the following 
values of parameters are chosen in the case of layered model, 
u^=0.55, 1.0, 0.55, 1.0, a = 0.55, h^ = 0.09 

and hg *= 0.1, 13 = 10.0. ' 

The vertical concentration profiles of pollutant 
are depicted in figs. (6, 1-6. 4). It is found that the 
concentration of pollutant decreases as downwind distance 
increases in both point and line source cases. From these 
figs., it is also noted that the concentration of air 
pollutant decreases due to deposition on the ground as well 
with chemical reaction parameter. Similar results have 
also been obtained in the case of three 1 ay effidt model (graphs 
are not shown) , 
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The effect of stepwise variation in the wind velocity 
on the concentration distribution of pollutant is shown in 
figs. (e.S-d.e) for h^^O.OO v^= 2.0/ 0j^=lO.O using 
following set of parameters. 


(a) One layer 

u^=l,0/ k:^^=i.o 

U2=1.0. 2=1.0 

(c) throe layer 

U3_=0.55/ U2=0.75, u^=1.0 

h^=0.1 


(b) two 1 ayer 

u^=0.55/ u^^l.O 


hg=0.l 


K^l=1.0/ 


K 




It is found that as the number of layers increases the 

concentration near the ground decreaffis. 

Thu effect of stepwise variation on diffusion coeffi- 
cient in z~dircction are also shown in figs. 6. 7-6.8, for 
h^=0.09/ v^= 2 . 0 /Pi=lOkO,usingthe following set of parameters 


(a) One layer (b) two layer 


u^=l-0 

K^,= 1.0 


Uj^=1.0/ u^^l.O 


h =0.1 

g 


(c) three layer 

u^=l,0/ U 2 -I.O/ 
h^=0.1/ 


tr n cc K =1 0 K 1=0,55/ K„,=0,60 

u =1,0/ 


K^3=1.0 


In this case the concentration of pollutant increases 


«T,mVirtr of layers increases , 


in the lower region as 
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CHAPTER VII 


STEADY STATE DISPERSION OF A REACTIVE AIR POLLUTANT IN 

PRESENCE OF GREEN BELT 

7, 1 INTRODUCTION 

As pointed out earlier the dispersion of air pollutants 
is governed by the processes of molecular diffusion and 
convection. In the atmosphere, dispersal process depends 
upon the type and number of sources, stack heights, various 
meteorogical factors (such as wind, temperature inversion, 
rainout/washout) and topography of the terrain (Pasquill, 

196 2; Seinfeld, 1975), Several studies have been conducted 
to understand the process of pollutant dispersion by 
including some of the above mentioned factors (Smith, 19 57; 
Pasquill, 1962; Hoffert, 1972; Slinn, 1974; Scriven and 
Fisher, 1975; Lamb and Seinfeld, 1973; Calder, 1977; ErmaX, 
1977; Alam and Seinfeld, 1981; Reda and Carmichael, 198 2, 

Llewelyn, 198 3/ Karamchandani and Peters, 1983). in 

particular, the redistribution of a gas plume caused by 
reversible washout has been investigated by Slinn (1974). 

The modelling of atmospheric pollution by nitrogen and 
sulfur dioxide; has been studied recently by Alam and 
Seinfeld (1981) and Reda and Carmichael (1982). Tbe deve- 
lopment of a second generation mathematical model 
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air pollution has been presented by Gregory et*al, (198 2) . 

The effect of foggy environment on reversible absorption 
of a pollutant from an area source has recently been 
studied by Shukla et-al. (198 2). 

The transport of gases and particiolate matter within 
plant and vegetable canopies has also been studied.Petit et«al. 

(19 76) presented results concerning characteristics of air 
flow within and above a forest by calculating SO2 fluxes 

at the top of canopies. Bache (1979a, 1979b) used modified 

/ 

form of diffusion equation to study particulate transport 
within and above the foliage canopy, Slinn (198 2) gave 
a theoretical framework to predict particle deposition due 
to vegetation by considering a variable wind velocity 
profile. A review of atmospheric deposition and plant 
assimilation of gases and particles has been presented by 
Smith (1957) and Hosker et al. (198 2) wherein a mathematical 
model for aerosol depletion and deposition on forests 
employs a modified form of convective diffusion equation. 

This model develops the interaction between forest 
structure and open field by considering forest aerodynamics 
and aerosol characteristics (Wiman# 1965)). 

From the above studies it may be speculated that if 
a suitable green belt in the form of forest is provided 
around tho pollutant source, some distance away from source 
but Close to tho place (e.g. a habitat, a historical monui<.ent) 
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to be protected, then it is possible to protect the endan- 
gered area under consideration, m this direction Kapoor 
and Gupta (1984) studied the attenuation of an inert pollu- 
tant by a green belt under steady state conditions. In 
this chapter, we study the dispersion of a reactive air 
pollutant from point and line sources. We focus on dry 
deposition on the ground when the wind velocity and diffusion 
coefficients are functions of elevation. The exact solution 
of the diffusion equation with reaction term is obtained by 
dividing the inversion layer into four parts. Wind velocity 
and diffusion coefficients are taken as step functions with 
smaller values in the lower region than in the upper region. 
The effect of green belt on the reduction of concentration 
of pollutant due to a removal mechanism is, therefore, 
discussed. 

7. 2 mathematical formulation and solutions 

(i) DISPERSION FROM POINT SOURCE 

Consider the dispersion of a reactive air pollutant from 
a point source as discussed in chapter VI. Assume that a 
green belt (a suitable tree plantation capable of absorbing 
pollutant) of thickness d is located at a distance from 

the source along the wind direction and extending in- y 
direction. The physical situation is illustrated in fig. 

7, 1 where the inversion region is divided into four regions, 
the green belt being in region III. The partial differential 
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cquatiore governing the concentration distribution of the 
air pollutant in different regions are written as follows: 

Region I (0 ^ 2 < h^) 


9C. 


u, 


3x 


K. 


a c, 


+ K , 

^2 yl 


3 C 


ay 


1 

i -kC^ 


(7.1) 


With boundary conditions 


(1) 

Cf = 0 

at 

X = 

0 

(ii) 

0 

H- 

ll 

0 

as 

y -* 

T 00 

(iii) 

ac, 

V n r— ~ 

Ki 3 2 

« V 

dl^i 

at 2 = 


3 a b 

(i'') -51; “ Ts ' b 


C7.2) 
(7. 3) 

0 (7.4) 

= at 2 = hg (7.5) 


Region li (h^ < 2 < H) 


u 


2 ax 


K 


z2 


3^C 

a 2 


2 + K 
2 Y2 


ay 


2-’^h 


and bounduiT conditions are 


(i) 

(ii) C 


6(y)6(2-hs^ ^ 

3 


0 


as y z ^ 


ac^ 

(Ixi) — T ^0 at 2 H 

3 ^ 


9 


a ^ 2 tr ^ 
(iv) 2 "Tz “ 9^ 


® K„ i ~n « ^ ^2 


C, at 2 = h 


(7.6) 

(7.7) 

(7.8) 

(7.9) 

( 7 . 10 ) 
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Region III (0<z<h) x>x 


u 

""3 3x 


K 


0 

3 C. 


^ 3^C 

2 3 . _ 2 ^3 -(k+X)c. 

yj ay2 , 


3z 


boundary conditions are 


(i) C, 


C T at X 


X-, 


(la) ^ as ^ 


3 C. 


3C, 


^-^3 "tI ' ' = 3 = C at z=h 


■g 


ac. 

Z3 3z ^2 "3 


" 3 

(i’v) iy, -, = v^, . C„ at 2 = 0 


Region IV (b < 2 < H) x > x. 


3C 


2 

3X, 


3^C, 


'4 -'4 A 

^4 ax'" ^z4 TTa ^y4 Tn 


- kc 


3z 


ay 


4 


with boundary conditions 


at X = x^ 


c , =0 as y ^ + 


00 


3C^ 3C- 

^z4 "Tf = ^23 TI ' V ^3 at z = 


1^4 
a 2 


= 0 at z = H 


(7.11) 

(7.12) 

(7.13) 

(7.14) 

(7.15) 

(7.16) 

(7.17) 

(7. 18) 

( 7. 19) 

(7. 20) 


whore C^, are concentrations of pollutant in 

different regions, k is rate of chemical reaction and X 
is depletion rate of pollutant due to green belt. 
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It is convenient to cast the problem in dimensionless 
form. To do so we use following dimensionless quantities 


- _ c _ z r, = ^i._ ~ _ Y 

^ ^ - K' u. - — 

max max 


H 


K 


K . 


Zi 


ai K, 


' "^di ■■ ^ ^ ' Ma K, 

zmax 5 

h 


'zmox 


_ h __ 


lad - _ ;d£ 

d.2 “ K„ 

max 2 max 

2 


s H 


. C. = 


(7.?1) 


u H ' 

- 4;- -'- C. (i=1^2,3,4) . 
W 1 


Equations (7.1—7,20) can be v/ritten in the 


dimensionless form_as 
Region I (0 < i < hJ 


3 C- 


X 


e, —A + T 

ay 


5^1 

^ 9 z ^ 


Cl = 0 


C3_ = 0 


at X = 0 


as y ± " 


Vi 


(7.22) 

(7. 23) 
(7.24) 


3 C, 


2 


= at z = 0 


(7. 25) 


K 


zi 


ac. 


3 z 


= K 


z2 


3C, 
3 2 


' C^- ^2 at z = h^ 


(7. 26 ) 


Region II (hg < z < 1) 


lls 

3x 


3 *7o 

|3 + y 


Or 

3'^ 


2 "2 ' '2 o g 

^ ay 3z 


2 - a c 

2 2''2 


(7. 27) 
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c, = r 6 (y) 6{z-h ) at X = 0 

^ u “ 

(7,28) 

as y^+oo 

(7. 29) 

8 C„ 

=0 at z = 1 

8 z 

(7. 30) 

ac _ 3 c. 

(7.31) 

Region 111 (0 j< z < h ) x > x. 

y X 


ac a^c a 

— J = 0 — ^ + r ~ (a +x )c- 

ax 3 ay^ 3 8 z2 3 3 3 

(7.32) 

at X - 

(7.33) 

C„ = 0 as y-ttoo 

3 "" 

(7.34) 

bS 

= N C at z = 0 

0 z ^ 

(7.35) 

3 C„ _ 9C _ - _ _ 

~ ^z4 ^3” ^4 ^ ~ 

^'•^az 8 z 

(7.36) 

Region IV (h„ < z < 1) x > x. 

g - - X 


DSi sy^ “952 

(7.37) 

X = X^ 

(7. 38) 

C^ = 0 as y -* ± “ 

(7. 39) 
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15 

3 z 


0 at z = 1 


K 


z4 


3 C, 


3 z 


3 C, 


K 


/ C.= c at 2 = h 


" '"4 -g 


where 


K 


0. = - 


K U. 
zmax 1 


K . 
Zl 


u. 

1 


N 


= N = ^ a - a - “ 

K^2 ^3 


^ “ r — Ah ^ 

j - ^ A ~ 

3 zmax 


(7.40) 
(7. 41) 


The solutions (x,y/Z), C^Cx/y^z) of equations 
(7.22) and (7*27) subject to boundary conditions (7, 23'*7. 26) 
and (7. 28-7.31) have been obtained in Chapter VI, 


Similarly, the solutiors of equations (7,32) and (7.37) 

under the condition , and subject to boundary 

^3 4 

conditions (7. 33-7.36) and (7.38-7,41) are obtained as 
follows (dropping bars for convenience) : 


00 (x-x^) 

C^(x,y,z) = - — E e 

3 n=l 


R« (z) (7.42) 

n In 






■n5'^ (X-Xj 


E e 


n 


N/li^ n=l 


l.R^G^^U) (7.43) 
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where. 


/ U 2 C^(x/ 2 )G^j^( 2 ) dz u^C 2 (x,z)G 2 j^(z) dz 

p I s P„ , ■■'■ — 

1 

^ G^^(z)dz + / G|^(z)dz 

o hg 


G.^(z) = 


9 

sin al.z + cos aJ^z) 

^11 li_ 

sin a^Li a‘ihg) 


(7, 44) 


(tan al ^ sin al „2 + cos 

Q (.^) = 

2n (tan a ^2 sin V ^iaV 


= 0,Cx-Xi) + 3iXj 


^ii= -Vh-^^’'3 ^12° 




and 31^® eigen values of following equation 


(N cos " ^11 ^il 


"23 N- 1 V, ^ 

(■^1 sin a^ 2 . ^11 “g 


tan a![ ^ cos h^- a]^^ sin a ^2 ^gL (7.45) 

(tan a 5_2 ^12 ^g ^^2 V 

When the source lies in the second region, the 
concentration distributions C 3 _(x,z), C^Cx^z) involved 
in equation (7,44) are given by equations (6.39) and 

(6.40) in Chapter VI. 
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Similarly when source lies in the first region# 
corresponding C^Cx^z) are given by equations 

(6,41) and (6,42) in Chapter VI, 


(ii) DISPERSION PROM A LINE SOURCE 


In this case, the concentration distributions 
C^(x/Z)jf C^Cx, z) corresponding to region I and II have 
been obtained in Chapter VI and are given by equations 
(6,39) and (6.40). For the region III, IV, the concentration 
distributions are obtained by integrating equations 

(7. 42“7. 43) between -oo to <» with respect to y and 
written as follows. 


*1 


C (x, 2 ) = S G,„(z> 


n=l 

00 

C (x,z) = S e 
n=l 


• 6' ^(x-xj 

R1 G^( 


n 2 n 


z) 


(7. 46) 


(7.47) 


where R^ is defined by equation (7,44) . 

7.3 RESULTS AND DISCUSSION 

To study the depletion of pollutant due to green- 
belt, the following valties of parameters have been chosen 
for computation: u^ = 0.55, U 2 = 1.0, U 2 = 0,5, u^- 1,0, H— 200m 
K^l= 0.55, 1.0, K ^3 = 0.5, = 1.0, = <3^ « 03 = 

p = 10, a = 0 , 11 , X = 4,0 and h = .05. By taking various 
heights of source^the vertical concentration distributions 
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of pollutant havebeen computed at the end of green belt 
and depicted in figs, 7, 2-7*4 It is noted that the 
concentration of pollutant decreases in the presence of 
green bolt. 

If we compare figs. 7.3 and 7.4 it is found that 
the depletion in the concentration of pollutant due to 
green belt is greater if the height of source is less than 
the height of green belt. 
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CHAPTER VIII 


UNSTEADY STATE DISPERSION OP A REACTIVE AIR POLLUTANT ; 

EFFECT OF GREEN BELT 

3*1 INTRODUCTION 

The effect of removal on the steady state dispersion 
of a reactive air pollutant due to presence of green belt has 
been discussed in Chapter VII » In this chapter, the same 
problem is investigated for unsteady state dispersion of a 
reactive pollutant in the atmosphere emitted from a point 
source# It is assumed that the green belt lies away from 
the source in the wind direction and is acting as a sink so 
that the pollutant can be removed by chemical reaction, 
deposition absorption^ etc# In the region of green belt, we 
consider the unsteady state diffusion equation with a removal 
term proportional to the concentration of pollutant such that 
the magnitude of removal rate coefficient depend on the number 
and size of plants in the green belt (silnn#1982/ Wiman and 
Agren (1985))# However, in the region without green belt; the 
usual unsteady state diffusion equation with chemical reaction 
i s taken • 

8.2 MATHEMATICAL FORMULATION AND SOLUTIONS 

Consider the unsteady state dispersion of a reactive 
pollutant from a time dependent point source in presence of gree 



123 


belt under the inversion condition# Assume that the green bel' 
of thickness d is located at a distance from the source in 
the wind direction* The region under consideration 0< 3c<x^+d 
is divided into two partS/ the green belt is assiimed to exist 
in the second region# It is assumed that the wind speed is 
sufficient large in both the regions so that diffusion in 
X“direction can be neglected in comparison to advection# In 
writing the mathematical model it is considered that the green 
belt is equivalent to point sinks distributed all over the 
second region# ThuS/ the partial differential equations 
governing the concentration of pollutant in the two regions car 
be written as follows : 


Region I x ^ 


aa ac. a^a 

* -1 anr = V ^ 

with boundary conditions 

= O at t = 0 

C, = 6(y) a(z-h 3 ) 

1 

= 0 as y •* ± 
ac. 

—a. = 0 at z = H 
az 



at X = 0 


ac. 

K f r— 
zl dZ 


(S#l) 


< 18.2 ) 

(S#3) 

(3.4) 

(8.5) 


vai at z = 0 


C3.6) 
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Region II 

X > 







__ a 

at- 

ac, 

+ U^ "T-" 

a 9x 

a _ ^^^a 

^va 2 

+ 

K 

za 

a^c 

— 2 

(8.7) 




ay 




az 


and 

boundary conditions 

are 





^a 

= 0 

at 

t 


0 


(3.3) 


^a 

= 

at 

X 




( 8 *9 ) 


^a 

= 0 

as 

y 

-> 

+ ' 

DO 

(8.10) 


3^a 







(3.11) 


az”" 

= 0 

at 

z 


H 










(8.12) 


^za 

Q ^ 

’'32 

at 

Z ; 

= 0 


where are concentrations of pollutant in each region, 

are mean wind velocities, ^za diffusion 

coefficients in y- and z-directions respectively* v^ are 

dry deposition velocity on the ground in the different regions* 
k is chemical reaction rate, \ is the deplection rate of pollutan 
due to green belt* 

As in Chapter IV, the following forms of W(t) would be 
considered in the subsequent analysis : 

(i) Flux is instantaneous 

w(t) = 6(t) 

(ii) Flux is constant 

W(t) = (constant) 


(3*13) 
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(iii) Flux is step function type 
W(t)=:W„ 0<t<t 

C — o 

= 0 t > t^ 

The concentration distribution of pollutant for region I is 
(i) For instantaneous case 


C. (x,y,z,t) = “2 P(x,y,z) 6(t - ^ ) 

1 

(ii) For constant case 


(3.14) 


C^(x,y, 2 ,t) = P(x,y,z) H(t ~ 

1 1 

(iii) For step function type flux 


(8.15) 


rV 

C. (x/y/2,t) = -- P(x,y,z) [H(t - ^) - H(t - ~ 

i U^ Uj, Uj 


where 


t ) H(t-t )] 
o o 


(3.16) 


exp(jJ— -) 
— iL. X 


P (x,y^z) 


/ K n=i ^n 

\j^ ^ naiiiil!ii Ill 

^1 


E cos Xj^(z-H) 


X cos X (h -H)e 
n s 


.(JS. flai Kh: 

“l '>1 “ 


n ^ 

P =5 / cos^ ?v (z“H)dz 
n n 

o 

and Xj^^s are solutions of following transcendental equation 
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V 




K 


ai 

2l 


(8.17) 


The concentration distribution of pollutant in second 
region is obtained as : 


(i) Flux is instantaneous 


W X. x-x. 

<=3 = 'GT - <G- + "GT^’ 

1 1 a 


(ii) Flux is constant 

W X. X-XL 

C = ^ M(x,y/ 2 ) H(t - (-- + rj — )) 
a "a 


(3.18) 


(3.19) 


(iii) Flux is step function type 

M(x,y, 2 ) CH(t-(i + i:j^))-.H(t-t^-(^ t -^)) H(t-t^)] 


w 

C = -£ 
a 


“1 




( 8 . 20 ) 


where 


<e« 

(- rlrrrr^ °° co® Mjj(2"^^) 


M(x,y,2) = — ;z';', 

8(xJ n=l 


Q, 


n 


X exp 


K 


[ -lc(^ + - — <x-=<i > - 

^ a a a 


X r cos X^^Chg-H) 65^ [- ^ X^] 

m=:l m 1 


X / cos Xjjj(2-h) cos )UIjj(2-H) dz 

o 

H 

Q = / cos^ M„(2 -h) dz, 
o 

u ar© roots of following ecjuation 
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II tan U H = ^ 
n n K 


d2 


and 


za 


X. K 

0(x) - — Z-p — + (x- 3 ^) 


Using the following dimensionless quantities 


X = 


''z max ^ 


max 


1 

H 


j,y = j^,2=| ,U, 


u. 


u_ 


U, 


max 


' ^a = 


^ax 


K 


K 


2l 


zi ” K 


z max 




U 

max 1 




23 ” 


'^ai H 

> 


K 

2 max 


z max 

U H^c 


, ^2 

max a 

k 

k H 


~ K 


, t = 


K t 

z max 


z max 


K - 


X F/ 


W H 
o 


K 


z max 
h 


^ ^O = W K 


c z max 


rr,f^\ w(t) 
^ W(t:) = 


hg — f = Xj^H f , (i=l/2/ n=i/2/»«») (8*21) 


the concentration distribution of pollutant given by equations 
(3^4 - 3.16) can be transformed in dimensionless form 
(dropping bars for convenience) as/ 


Region I x < Xj_ 

(i) Flux is instantaneous 

w 

Ci{x,y#2,t) =~P(x/y,z) 6(t-‘J-) 


u. 


(ii) Flux is constant. 


C. (x,Y, 2 ,t) = P(x/y,z) H(t - ~) 
■^1 1 


(3#22) 


(8*23) 
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(iii) Flux is step function type 
P(x/y,z) ^ ^ _ 


(x,y,z,t) 


[ H(t-^) - H(t-t - ^^)H(t“t )] 


where 


exp(- 


(3*24) 


P(x,y,z) = 


z max 


/~ K , 

zl X 

z max 1 


“ 1 

S ^ \ (z-X) 


" p 

n=l 


cos Xn^hg-i) e 


.(JL x2)^ 


(8.25) 


P^ = / cos"^ Xjj(z”l) 
o 

and Xjj^s 3^® roots of following transcendental equation 


^n ^n “ 




Region II x > x^ 

(i) Flux is instantaneous 


X. x-x. 

U. + ^)) 


” IT t ■'■ "u" 

X la 


(8.26) 


(ii) Flux is constant 


M(x,y,z) H(t ~ + 


x, x-x. 

+ _J1)) 


(8.27) 
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(lii) Flux is step function type 
M(xiY/Z) _ 


^ [H(t - + - 5 ^)) - H(t-t„- ^)) H(t-t^ 

JL X at X 3 


( 8 . 23 ) 


where 


°° COS JU,„ (2"i ) 

M(x.y,z) s 


n 


'/471 / 3 (x) n=l 




n 


K 


X exp [-k(^ + -^) - ^ (x-x^) - - 2 l ^j^(x-x^)] 


K 


X E ~ cos X„(h«-l) exp [- 


ms 

m=x m 


X / cos X^iz-l) cos iXj^(z-l) dz 


{ 3 , 29 ) 


Qjj == / cos*^ Mj^(z-l) dz 
o 


0(x) = 3j^Xj^ + ^g(x-x^) 


K 




1 “ K 


ii. 


z max 




K. 


/ ^ 


ya 


® max 


"'z max = 


u^^^ = maxCuwU ) 
max A a 


and ILL * s are roots of following equation 


'n 


'd 2 




n K 


za 
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STEADY STATE CASE 


In such a case t -* <», the corresponding expressions 
for concentration distributions in the two regions can be obtaine 
from equations (S#23) and (3*27) as follows (in the case 
of continuous point source) : 


_ X p(x,y,z) 
= M(x,y,z) 


where P(x/y/z) and M(x/y/z) are defined in previous section 
[see equations (8»25), (8ii29)] • It is remarkable that these 
expressions can also be obtained from equations (6*39 •• 6*40) 
in Chapter VT assuming u^ * ^2^ ^zl “ ^z2' ^yi “ ^y2 
taking U 3 = u^ = u^, K ^3 = K ^3 = 

equations (7*42 - 7-43) in Chapter VII* 


K in 
ya 


8.3 RESULTS AND DISCUSSION 

For unsteady case# the concentration distributions of 
pollutant for different values of t# 

e, = = 10 * 0 . a = 0 * 11 # = 5 . 0 # are computed and shown in 

1 a 

figs. (8.1#3.2#8.3)* It is seen fron these graphs that the 
concentration of pollutant decreases in the region of green belt 
in all the three cases. The vertical concentration profile of 
pollutant is depicted in fig. (3.4) for K = 6 . 0 # 10.0# u^^ = 1.0, 
~ ^zl “ ~ ^ * 10.0# = 10.< 

and x^ = (end of the green belt). It is observed t.hat the 
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concentration of pollutant decreases as X increases* 

The analysis presented in this chapter shows the 
importance of green belt in reducing the concentration of 
air pollutant emitted from time dependent sources such as 
accidental leakage of gases etc. and there by protecting 
region under consideration* 
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CHAPTER IX 


ATMOSPHERIC DISPERSION MODEL WITH INTEGRAL 
TERM FOR REMOVAL MECHANISM 

9 . 1 INTRODUCTION 

The Study of atmospheric dispersion of air pollutants 
from point source/ line and area has received a great deal 
of attention during last few decades (Smith/ 1957? Pasquill/ 
1962; Gifford -and Hanna/ 1971/ Heines and Peters / 19733/ 19731)./ 
Lebedeff and Hameed/ 1975; Ragland/ 1973; Ermak/ 1977; 

Dobbins, 19 79; Karamchandani and Peters, 1983). Effects of 
removal process on the dispersal of air pollutants have 
also been investigated (SQriven and Fisher, 19753, 19 75b, 

Alam and Seinfeld, 1981, Fisher, 198 2, Llewelyn, 198 3). When air 
pollutants are removed by rainout/washout in the environment 
the dispersal process becomes very complex and very little 
attention has been given to understand this mechanism 
(Hales, 1972; Slinn, 1974; Alam and Seinfeld, 1981). 

In view of this, in this chapter a growth equation is 
suggested to study the removal of air pollutant by rain or 
fog droplets. In such a case the -unsteady state partial 
differential equation for dispersion of air pollutant is 
modified into the form of an integro- differential equation. 


■^Published in Applied Mathematical Modelling Vo. 11, 1987 
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To explain the effect of removal process by this model, 
we have studied the dispersion of air pollutant from an 
elevated time dependent point source which may be applicable 
during mansoon season in India. 


9.2 MATHimATICAL FORMULATION 

Consider the dispersion of air pollutant from an 
elevated time- dependent point source in a stable environment 
such that the wind velocity may be assumed to be laniform. 

In such a case, the three dimensional transport equation 
in cartesian coordinates (x,y/Z) for air pollutant concen- 
tration C is written in the following form 


+ u 

9 1 9x 


D(^ + ^ ^ - 

3^2 ay^ 


kC 


C9.1) 


where U is mean wind velocity, D is diffusion coefficient 
assumed to be constant and k is the removal rate. 

vmen the fog or rain droplets are present in the 

environment, the air pollutant gets absorbed in 

, . TM-nr'pqs is not instantaneous, 

droplet x^hase and this removal pr 

that tha rate of decrease of pollutant (i.e. 5 ^) 
is proportional to its concentration in the droplet 

^ ^rYii;:,-hion (9.1) 

phase, the transport equ 
modified as 

s •" g -A • • S’ ■ 
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Since the rate of change of the concentration of air 

3 C 

pollutant in the droplet phase (i.e. may increase 

with the increase of concentration of pollutant C and 
decrease due to rainout/ washout/ the governing equation 
for can be written as. 


3C 

= a C - ac 
at o r 


(9. 3) 


where constant rates of absorption and removal. 

Now integrating equation (9.3) with the condition = 0 
at t = 0/ we get 


C^Cx/Y/Z/t) = a / e”“^^“'^^G(x,y/Z/.T) 


dl 


(9.4) 


Equation (9.2) then becomes 


+ TT iS - 

at 3x “ 



ais . / e-“‘^-’^'c(x,Y,.,T)dT 

ay^ 3z^ o 

(9,5) 


where the integral term represents the removal of air 
pollutant by the droplet phase and 

It may be noted here that equations (9,1) and 

(9,5) can be combined as follows: 


2 2 2 t 

^ ^ ^ "" D(^ + + L^)- / G(t-T)C(x/y/Z,T)d!r (9.6) 


3 t 


3x‘ 


3y‘ 


3 z ' 


where for G(t-'T) = k6(t-T)/ 6(,) is Dirac delta function, 

cec 

we get equation (9.5) . 


-tt(t-T) 

the equation (9.6) reduces to (9.1) and for G = a^e ^ 
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Thus^ equation (9.6) may be interpreted as a general 
transport diffusion equation with an integral term represen- 
ting the removal mechanism. Equation (9,1) may be designated 
as diffusion equation with instantaneous removal term while 
equation (9,5) may represent the case with delayed (slow) 
removal process. 

It is noted here that under steady state condition 

a k 

equation (9,3) gives — C. In such case^ 

equations (9,1) and (9,2) would become identical provided 
a k 

k = — s= ~ and two removal mechanisns would be same under 

this condition. However^ under xonsteady condition these 
removal mechanisms are different. 

In the following, therefore, we study the dispersion 
of air pollutant with delayed and instantaneous removal pro- 
cesses under the condition when k is greater, less or 
equal to 

The initial and boundary conditions for (9.1) and 
(9.5), taking the source as the origin of coordinate system, 
are written as 

C(s,t) =0 at t = 0 for all s = x^+ y^+ z^> 0 (9.7) 

C(s,t) =0 as s -» CO for t > 0 (9.8) 

“471 s^D ~ = W(t) as s -» 0 for t > 0. (9.9) 

8 s “ 

The last boundary condition implies that the point source 



136 


has prescribed time dependent flux. In this chapter^ the 
following forms of W(t) are considered, 

(i) W(t) -■ W^6(t) 

(ii) W(t) = 

c 

(xii) W(t) = W 0 < t < t (9.10) 

W wmm Q 

= 0 t > t 

o 

where 6(«) is Dirac delta function. 


Define the following dimensionless quantities 

..2 


t = 


u 


t, U = 


U- a = a 


n _ D „ - u* - U'^ - _ U'^ 

a - _ a . s = 5-s, X = ^ x,y - 5-y 


(9.11) 


U 


liTz k = JSS- 

Z/ Jv — ^ f 


D 


u^c-2 


C = 


u*W 


D 


C. 


where u'^ is friction velocity. 

Using above dimensionless quantities equations 

(9,5~9.10) can be written in dimensionless form (dropping 
bars for convenience) as 


^ + u 

3C _ 

a^c 

? 9 

3 C 3'^C 

at 

3x 


9 9 

dY 3z 

C(s,t) 

= 0 

at 

t = 0, s > 0 

C(s,t) 

= 0 

as 

S *+ CO t ^ 0 


/ e ®^^“^^C(x,y,z,T) dT 
o 

(9.12) 

(9.13) 


(9.14) 
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.47,s 2 as s 

as 

c 


0 t > 0 


(9.15) 


(i) Slil = Q„6(t) 


Ui, 

c 

(iii) =1 0 ^ < "^Q 


(9.16) 


0 t > t 


where 


o 


W D 

U \ 


9.3 METHOD OP SOLUTION 

Talcrng Laplace transform of (9.12) and using initial 
condition (9.13), we get 


„ 3C _ 3^0 a^^ 3^ ^ 

^ 3x^ ay^ 32 *^ 


(9.17) 


a. 


where C is Laplace transform of m = p + and p is 

Laplace variable. 


The boundary conditions become 
C(s,p) - 0 s 


-471 s2 ^ s ^ 0 

as • 


(9.18) 


Assximing the solution of equation (9.17) of the 


form, 
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C = exp (~) f(s,p) 
and using (9, 18) , we get 


c - 

47tSW 

c 




+ m) ] 


(9.19) 


Taking inverse Laplace transform of equation (9.19)/ 
the solution of (9.12) satisfying the initial and boundary 
conditions can be written as follows 


C(x/y/Z,t) = 


1 

47ts 


(— 1 _) 


T+ioo 

/ 

y_ioo 


w(£i 


w. 


exp [pt + 


2 


where 

(a/a*) = •” (~ + a) + {(^l - a) 2-4a^}V2 and r is a 

real positive number such that all the singxalarities of the 
integrand lie on the left handside of line Re(p) = ^ in 
the Bromwich contour (Carslaw and Jaeger/ 19 59) . The 
integral in equation (9.20) can be evaluated and the concen- 
tration along the central line for each case can be found 
as follows. 


(i) when the flux is instantaneous/ i.e. W(t) = W^6(t) 

Q P(x) «• 

C(x/0/0/t) = — / M(x,U/t)du + / M(x/U/t)du (2.21) 

471 '^x a a' 


where P(x) = exp (Ux/2), 
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M(x# u, t) 



Similarly the solution of equation (9.1) is obtained and 
the concentration along the central line in the dimension- 
less form is written as 


C^(x/0,0,t) = 


(4T!;t) 


exp [ 


(9.22) 


when U = 0 and k = 0/ it is reduced to same t' form as the 
one obtained by Carslaw and Jaegar (1941). 


(ii) When the flux is constant^ i.e. W(t) = 

2 Oi* 1/ 2 

C(x,0^0 ,t) = {exp [-(|- + X -jj 

- ^ ~ M(x,U/t) du 

a 

1 CO -1 

- “ / ^ M(x,u,t)dUl 


(9.23) 


From the above expression it is observed that the concen- 
tration of pollutant as t ^ reaches to steady state/ 
given by the following expression 

c(x,0,0,t, 


Similarly solving equation (9.1) in this case the concen- 
tration dustribution along the central line in the dimension- 
less form, is given by 

C^(x/0/0/t) = C-(f“ + k)^/2^ ] 

«. 1 1/2 

- —f — e3q3(-ut) sin(u- a^) xdu} (9,24) 

^1 



140 


where 


a . = 


U‘ 


1 4 


(ixi) When the flux is step fiinction i.e, W(t)=VJ 0 < t < t 


=0 


t > t 


o 


P ( x) r 2 

C(x,0,0,t) = |g^£exp[ -(| + ^) x] (l-HCt-t^)) 


. a ut 

^ - [l-e H (t-t^)]M (x, u,t) du 

a 


ut 


1 

f “ [l-e H(t“t )] M(x/U/t) duj (9,25) 


7t U 

a* 

Where H(t-t^) is defined by (Carslaw and Jaegar, 1941) 


H(t-t^) = 0 


t ^ -^o 


t +€ 
o 




=1 t > 

o 

In this case, for instantaneous removal the concentration 
distribution along the central line, in the diemsnsic.nless 
form, is given by 

C^(x,0,0,t) = {exp [-(^ + k)^/2x] X 

[l-H(t-t^)] - ^ ~ e'^^sin [(u-a^) ^^^x] x 

^1 

ut 

[l-e °H(t-t )] du}. 


(9. 26) 
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where 


a . = 


U' 


(iii) When the flux is step function i,e. 


W(t)=VJ^ 0 < t < t^ 


=0 


t > t. 


c(x,0,0,t) _ l^fexpf -(| + ~) ^ X] (i-H(t-tJ) 


1 1 


ut. 


~ u H (t-t^)]M 


1 oo . Ut 

5 ^/ - [1-e °H(t*-t^)]MU,u,t)du3 (9.25) 


Where H(t-t^) is defined by (Carslaw and Jaegar, 1941) 


H(t-t^) = 0 


_ t 


•fc irG 
O 


^ -S 




= 1 


t > t 


o* 


In this case, for instantaneous removal the concentration 
distribution along the central line, in the (3i©msnsic-nless 
form, is given by 

2 

Ci(x,0,0,t) = (exp [-(f + k)V2x] X 


^ r j e""\in [(iwa ) x 


Ut 


[l-e °H(t-tJ] du}. 


(9.26) 
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Now to compare the case of delayed removal with 
instantaneous removal we consider the following cases 

a. a a 

(i) k > — (ii) k = ^ (iii) k < 

(X 

1 

For k > ~ , it is observed that for all the three types of 

sources, the concentration of pollutant for the case of 

delayed removal is greater than in the case of instantaneous 

removal (see figs, 9.1, 9,4, 9,7), In the case of constant 

1 

flux, it is noted that even if k _< , the concentration 

of the pollutant for the case of delayed removal is always 

greater than instantaneous removal case (see figs. (9. 4-9. 6)), 

a 

1 

However, for other two typesof sources when k < the 
concentration in the delayed removal case can be less, equal, 
or greater than the instantaneous removal case and this 
behaviour depends upon the downwind distance from the source, 
time elapsed, etc. (see figs. 9. 2-9, 3, 9, 8-9. 9) , The 
results can be visualized by comparing the removal 
mechanism in equations (9.1), (9,5) and noting that the 
integrand in the integral term in equation (9.5) involves 
a negative ej<ponential which decreases as t or ft increases. 

To observe the effect of wind velocity on the concen- 
tration of air pollutant at different locations and times, 
see figs, 9,10-9.12. It is noted from fig. 9,10 that, for 
instantaneous flux, the concentration of pollutant at a 
fixed instant increases at different locatiore and the point 
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of maximum in the concentration distance profile moves 
away from the source as wind velocity increases. For 
constant flux, (see fig, 9.11) the concentration increases 
as wind velocity increases at a particular time and 
location. However, when the flux is given by step function, 
(see fig, 9,12) the concentration increases as wind 
velocity increases for fixed t < t^, but when t > t^, 

the increase or decrease in concentration with respect to 
wind velocity depends upon downwind distance from the 
source. 












uoi|oj)u«3uo3 ssaiuoisuauiiQ 


Otmcnalonlesa downwind distonce 0Jm«n»(onte»8 downwind diatancc 

FIG NO G.7 FLUX IS STEP FUNCTION TYPE AT THE SOURCE. STEP FUNCTION TYPE AT THE SOURCE 







Dimensionlest downwind disionc* Dimcn#ionl8ss downwind dtstanc* 

rie NO 9.12 FLUX IS STEP FUNCTION TYPE AT THE SOURCE NO 9 11 FLUX IS CONSTANT AT THE SOURCE 
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